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Abstract 


Wave  packet  propagation  methods  are  used  to  compute  scattering  Wigner  Distribution 
Functions  for  the  square  well  potential,  square  barrier  potential,  the  H  +  H2  and  OH  +  CO 
molecular  collisions,  and  several  M  +  Ng  atomic  collisions.  The  scattering  Wigner 
Distribution  Functions  are  used  to  interpret  how  probabilities  flow  among  various  potential 
energy  surfaces  as  a  function  of  time  during  a  collision.  Positive  values  of  the  scattering 
Wigner  Distribution  Function  correspond  to  the  addition  of  probability  to  scatter  into 
the  state  that  corresponds  to  the  asymptotic  limit  of  the  potential  energy  surface  of 
interest.  Negative  values  correspond  to  the  loss  of  probability  to  scatter  into  the  state  that 
corresponds  to  the  asymptotic  limit  of  the  surface  of  interest,  and  zero  values  correspond 
to  probability  associated  with  the  wave  packet  that  is  still  in  the  interaction  region.  The 
loss  of  probability  on  one  surface  corresponds  to  the  addition  of  probability  on  another 
surface  at  different  times.  Bands  of  oscillating  peaks  and  valleys  that  form  in  the  structure 
of  the  Wigner  Distribution  Function  correspond  to  presence  of  secondary  transmission  or 
reflection  with  significant  probability.  The  square  well  frequencies  at  which  probability 
arrive  in  a  scattering  channel  corresponds  to  the  depth  and  width  of  the  well.  Scattering 
Wigner  Distribution  Functions  were  computed  for  the  following  combinations:  K  +  He, 
K  +  Ne,  K  +  Ar,  Rb  +  He,  Rb  +  Ne,  Rb  +  Ar,  Ci’  +  Ne,  and  Cs  Ar.  The  scattering  Wigner 
Distribution  Function  revealed  that  as  the  mass  of  the  noble  gas  increased,  a  significant 
proportion  of  probability  was  transferred  from  the  pump  state  to  the  lasing  state 
for  a  larger  number  of  total  angular  momentum  values.  Similarily,  as  the  mass  of  the  alkali 
metal  decreased,  there  was  a  reduced  transfer  of  probability  to  make  a  transition  from  the 
state  to  the  ^Pi  state.  The  reduced  probability  to  make  a  transition  from  the  ^Pi  to  the 

2  2  X  ^  2 

^P\_  manifolds  for  the  M  +  He  collisions  is  compensated  by  the  large  average  velocity  of 


IV 


He.  As  a  result,  the  K  +  He  system  exhibits  a  large  collisional  transition  rate  despite  the 
kinetic  bottleneck  caused  by  the  light  mass  of  He. 
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6. 1  Summary  of  the  M+Ng  collision.  The  states  1 3 , 3  >  and  belong  to  the  ^Pi 

2  2  2  2  ^ 

1 

manifold  in  the  asymptotic  regime  while  \\,\>  belongs  to  the  ^Pi  manifold. 
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Although  all  of  the  states  belonging  to  the  maniforld  are  at  the  same  energy 
level  in  the  asymptotic  regime,  they  were  artifically  separated.  Three  product 
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WIGNER  DISTRIBUTION  FUNCTIONS  AS  A  TOOL  FOR  STUDYING 
GAS  PHASE  ALKALI  METAL  PLUS  NOBLE  GAS  COLLISIONS 

I.  Introduction 

A  laser  has  three  properties.  The  first  property  is  that  there  must  be  an  active  medium 
that  interacts  with  photons.  Second,  a  population  inversion  of  the  active  medium  must  be 
achieved  in  order  to  achieve  gain.  The  population  inversion  occurs  when  more  constituents 
of  the  active  medium  are  at  a  higher  energy  level  relative  to  a  lower  energy  level.  Because 
a  population  inversion  doesn’t  occur  naturally  for  long  durations  of  time,  energy  must  be 
provided  to  the  laser  system  from  an  external  source.  The  process  of  exciting  constituents 
of  the  active  medium  into  a  higher  energy  level  is  known  as  pumping.  When  a  constituent 
is  excited  to  a  higher  energy  level,  the  constituent  will  spontaneously  emit  a  photon  and 
de-excite  to  a  lower  energy  level.  This  spontaneously  emitted  photon  will  travel  through 
the  active  medium  and  interact  with  the  other  constituents.  If  there  is  another  constituent 
at  the  same  energy  level  as  the  previous  constituent  that  spontaneously  emitted  the  photon, 
then  the  photon  can  stimulate  the  new  constituent  to  emit  a  photon.  This  stimulated  emitted 
photon  will  travel  in  the  same  direction  with  the  same  phase  as  the  spontaneously  emitted 
photon.  This  process  is  repeated  as  long  as  the  photons  aren’t  lost  to  the  laser  system.  This 
implies  that  the  third  property  that  a  laser  must  have  is  an  optical  feedback  system  in  order 
to  maintain  the  stimulated  emission  process.  Diode  Pumped  Alkali  Lasers  Diode  Pumped 

Alkali  Laser  (DEAL)  has  an  active  medium  that  consists  of  an  alkali  metal  generally  in  the 
gas  phase.  As  shown  in  Fig.  1.1,  the  alkali  atoms  are  in  their  ground  state,  the  i  energy 
level.  The  ground  state  alkali  atoms  are  excited  by  laser  light  emitted  from  diode  lasers 
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Figure  1.1:  Energy  Level  Diagram  for  an  arbitrary  DPAL  . 

(hence  the  name  DPAL)  into  the  manifold.  Through  non-radiative  processes  such  as 
collisions,  the  alkali  atoms  are  quickly  de-excited  to  the  ^P\_  energy  level.  Because  the 
energy  level  is  a  meta-stable  state  for  the  alkali  atoms,  a  build-up  of  alkali  atoms  can 
occur  on  the  from  the  ^Pi  manifold.  If  the  atoms  are  pumped  fast  enough  to  the  ^Pi 
manifold  and  are  quickly  de-excited  to  the  state,  then  a  population  inversion  between 
the  ^Pi  and  i  energy  levels  can  exist.  If  this  population  inversion  exists,  then  lasing 
can  be  achieved  between  the  and  ^  energy  levels  at  a  frequency  determined  by  the 
difference  in  energy  of  these  two  states. 

One  potential  problem  can  occur  when  the  non-radiative  processes  isn’t  quick  enough 
to  produce  a  build-up  of  atoms  on  the  ^Pi  energy  level  before  the  atoms  on  the 
spontaneously  decay  to  the  ^  energy  level.  In  order  to  maximize  the  output  power  of 
a  laser,  we  need  to  maximize  the  number  of  photons  that  underwent  stimulated  emission. 
This  implies  that  we  want  to  maximize  the  number  of  atoms  in  the  ^P^  level  relative  to  the 

1  level.  Therefore,  we  want  to  minimize  the  amount  of  time  that  it  takes  for  an  atom  in 

2 

the  ^Pi  energy  level  to  decay  to  ^P^  energy  level. 
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The  primary  non-radiative  process  that  de-excites  an  alkali  atom  from  the  energy 
level  to  the  energy  level  is  a  collision  with  a  buffer  gas.  Various  buffer  gasses  have 
been  experimentally  examined  as  potential  collision  partners;  however,  using  a  noble  gas 
atom  as  a  buffer  gas  is  of  interest  to  us  because  the  noble  gas  does  not  undergo  a  chemical 
reaction  with  the  alkali  atoms  which  depletes  the  number  of  atoms  available  as  an  active 

medium.  We  are  interested  in  the  rate  at  which  alkali  atoms  in  the  ^Pi  level  are  de-excited 

2 

to  the  energy  level.  This  rate  is  related  to  the  probability  that  an  alkali  atom  on  the 
degenerate  ^P^  energy  level  before  the  collision  will  end  up  on  the  energy  level  after 
the  collision.  For  this  reason  we  will  look  at  the  scattering  of  an  alkali  metal  (M)  with  a 
noble  gas  {Ng)  at  the  quantum  level  in  order  to  determine  which  noble  gas  maximizes  the 
probability  transfer  to  the  energy  level. 

Previous  studies  of  a  collision  involved  the  use  of  the  correlation  function  to  produce 
s-matrix  elements  via  the  Channel  Packet  Method  (CPM)  [8,  9,  11,  16,  17].  The  correlation 
function  is  a  time  signal  that  is  produced  when  the  wavefunction  of  an  evolving  reactant 
(initial  constituents  to  a  collision)  is  compared  with  the  wavefunction  of  a  stationary 
product  (final  constituents  of  a  collision).  The  s-matrix  elements,  which  contain  all  of 
the  information  of  a  scattering  event,  can  be  related  to  the  probability  that  a  reactant  with  a 
given  energy  ended  up  as  the  product.  The  problem  with  this  is  that  we  do  not  know  how 
the  probability  was  transferred  between  the  various  possible  states  that  a  reactant  could 
end  up  in  during  the  collision.  The  tool  that  we  will  use  to  accomplish  this  task  involves 
the  scattering  Wigner  Distribution  Function  (WDF).  The  scattering  WDF  is  the  Fourier 
transform  of  the  autocorrelation  of  the  correlation  function  with  respect  to  a  time  delay. 
The  result  of  this  Fourier  transform  is  a  quasi-probability  distribution  as  a  function  of  time 
and  energy. 

Although  the  Wigner  Distribution  Function  has  been  well  studied  in  the  past  [2,  7,14, 
15,  21],  the  application  of  the  WDF  as  a  tool  in  analyzing  a  collision  is  new.  Previous  work 
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involving  the  use  of  the  WDF  as  a  tool  to  analyzing  a  collision  examined  the  scattering 
of  a  reactant  in  a  finite  square  well  potential  [7].  The  results  from  this  work  can  be  used 
to  compare  the  accuracy  of  any  code  that  generates  a  scattering  WDF  for  a  square  well 
potential.  For  this  reason,  scattering  WDFs  were  produced  for  the  square  well  problem. 
Another  reason  for  producing  scattering  WDFs  for  the  square  well  problem  is  to  gain 
further  insight  into  the  meaning  of  the  scattering  WDF.  Lacy  was  interested  in  comparing 
the  scattering  WDF  with  other  signal  processing  techniques  such  as  a  spectrogram  and 
determining  what  features  of  the  WDF  can  be  related  back  to  the  parameters  of  the  square 
well  potential  or  the  transmission  coefficient  [7].  Because  our  work  focused  on  the  transfer 
of  probability  among  the  various  states  in  an  alkali  metal  plus  noble  gas  (M  +  Ng)  collision, 
we  will  use  the  simple  case  of  the  square  well  potential  to  interpret  the  meaning  behind  the 
values  that  the  scattering  WDF  and  how  the  structure  of  the  scattering  WDF  can  be  used  to 
determine  where  or  how  probability  was  transferred  during  the  collision. 

Once  the  scattering  WDF  for  the  square  well  are  produced  and  analyzed,  scattering 
WDF  are  computed  for  the  finite  square  barrier  potential  which  has  never  been  done  before. 
The  reason  for  this  is  that  the  potential  involved  in  the  M  +  Ng  collision  contains  features 
that  can  be  roughly  estimated  by  a  square  barrier,  scattering  WDFs  are  also  produced 
for  the  H  +  H2  and  OH  +  CO  collisions.  The  correlation  functions  used  to  produce  the 
scattering  WDFa  for  the  square  well  were  generated  by  code  that  I  developed;  however,  the 
correlation  functions  for  the  M  +  Ng  system  of  interest  were  developed  by  Loper  [9].  The 
H  +  H2  and  OH  +  CO  correlation  functions  were  also  computed  by  Weeks  [11],  and  the 
accuracy  of  the  code  in  producing  scattering  WDFs  from  correlation  functions  collected  by 
other  sources  can  be  compared  with  the  previously  produced  scattering  WDFs  for  the  one 
case  of  H  +  H2.  In  addition,  these  collisions  also  allow  us  to  confirm  our  understanding  of 
the  scattering  WDF  developed  from  the  square  well  potential. 
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Finally,  scattering  WDFs  are  computed  for  8  different  combinations  of  alkali  metal 
plus  noble  gas  collisions.  From  these  scattering  WDFs  we  hope  to  find  the  combination 
of  M  and  Ng  such  that  the  transfer  of  alkali  atoms  from  the  energy  level  to  the 
level  is  optimized.  Because  the  energy  level  is  degenerate  when  there  are  no  other 
collisional  partners  present,  there  will  be  several  states  that  the  probability  can  flow  to 
other  than  the  states  associated  with  the  energy  level.  Some  of  these  states  do  not 
allow  for  a  transfer  of  probability  directly  to  the  states  associated  with  the  energy 
level,  and  we  are  interested  to  know  if  the  probability  that  ends  up  on  these  states  in  the 
manifold  can  quickly  find  its  way  to  the  states  associated  with  the  energy  level 
or  if  there  is  only  a  one  time  transfer  of  probability  among  the  various  states.  Because 
the  strength  of  the  forces  that  allow  probability  to  transfer  among  these  states  changes 
as  a  function  of  total  angular  momentum,  we  will  produce  scattering  WDFs  at  various 
total  angular  momentum  values.  From  our  analysis  of  these  scattering  WDFs,  we  hope  to 
provide  conclusions  as  to  which  combination  of  alkali  metal  and  noble  gas  can  optimize 
the  lazing  performance  of  the  DPAL  system  that  other  sources  of  scattering  information 
such  as  the  evolving  wavefunction,  correlation  functions,  or  the  s-matrix  elements  can  not 
easily  provide. 
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II.  Scattering  Theory 


2.1  Scattering  Operator 

A  collision  consists  of  reactants  that  approach  each  other  and  interact  in  some  way. 
When  the  reactants  are  infinitely  far  apart  from  each  other  (the  asymptotic  region),  the 
individual  reactants  do  not  exert  a  force  on  each  other.  Instead  the  only  forces  present 
at  the  location  of  a  reactant  are  the  forces  that  make  up  the  reactant.  These  forces  can 
be  described  by  a  Hamiltonian,  and  the  sum  of  the  Hamiltonians  for  each  reactant  in 
this  regime  make  up  a  Hamiltonian  designated  as  the  asymptotic  Hamiltonian,  H^.  As 
the  reactants  approach  each  other,  the  reactants  exert  a  force  on  each  other.  This  force 
has  an  associated  potential  that  is  called  the  interaction  potential.  The  region  where  the 
interaction  potential  is  non-zero  is  called  in  the  interaction  region.  When  the  interaction 
potential  is  summed  with  asymptotic  Hamiltonian,  a  new  Hamiltonian,  H  is  created,  and 
this  new  Hamiltonian  is  called  the  full  Hamiltonian  since  it  includes  all  forces  involved 
in  a  collision.  After  the  reactants  interact  with  each  other,  the  states  that  result  from  the 
collision  are  called  product  states  [4].  If  lif/reactant  >  and  \iff product  >  are  the  initial  quantum 
states  of  the  reactants  and  products  respectively,  then  the  time  evolution  of  these  quantum 
states  in  the  asymptotic  regions  is:  reactant  >  and  U^\il/ product  >  where  the  time  evolution 

operator  in  the  asymptotic  region,  U^,  is 

=  (2.1) 

Now  let  =  0)  >  be  an  initial  quantum  state  subjugated  to  the  full  Hamiltonian  of  the 
scattering  event.  The  vector,  |i/r(t  =  0)  >,  is  called  the  Mpller  state.  The  time  evolution 
of  the  Mpller  state  is  =  0)  >  where  U  is  the  time  evolution  operator  of  the  full 

Hamiltonian  which  includes  the  interaction  potential.  The  Mpller  state  is  of  interest  for  the 
scattering  problem  because  the  Mpller  state  results  from  an  intermediate  quantum  state  that 
started  in  the  asymptotic  region  and  propagated  to  time  t  =  0  under  the  full  Hamiltonian. 
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The  reactant  state  is  also  created  from  this  intermediate  state  located  in  the  asymptotic 
region  except  the  intermediate  state  was  propagated  under  the  asymptotic  Hamiltonian.  In 
the  infinite  past,  the  quantum  states  =  0)  >  and  \i^reactant  >  must  be  indistinguishable 
(came  from  the  same  intermediate  state),  and  similarly,  in  the  infinite  future,  \i^(t  =  0)  > 
and  \iff product  >  must  be  the  same.  The  mathematical  statement  that  is  equivalent  to  this 
concept  is 


lim  \\U\lJ/{t  =  0)  >  reactant  >  II  =  0 

t—^—oo 

(2.2) 

\im\\U\il/(t  =  0)  >  product  >  II  =  0. 

t—^oo 

(2.3) 

This  condition  imposes  a  restriction  that  the  interaction  potential  tends  to  zero  faster  than 
All  interaction  potentials  used  in  this  study  meet  this  requirement.  The  Mpller  state, 
=  0)  >,  can  be  solved  for  in  Eq.  2.2: 

|(/^(?  =  0)  >=  lim  UW'^\l//reactant  >=  reactant  ^  (2.4) 

t—^—oo 

—  0)  >—  lirn  U  \4^ product  product  ^  •  (2.5) 

where  Q+  and  are  the  Mpller  In  and  Mpller  Out  operators  respectively  as  defined  by 
Eq.  2.4.  The  Mpller  operators  produce  a  Mpller  state  from  an  initial  state.  As  seen  by 
Eq.  2.4,  the  initial  reactant  (product)  is  propagated  backward  (foward)  into  the  infinite  past 
(future)  under  the  asymptotic  Hamiltonian  creating  an  intermediate  state.  This  intermediate 
state  is  then  propagated  forward  (backward)  in  time  under  the  full  Hamiltonian  until  t  =  0 
The  intermediate  state  created  by  the  Mpller  In  operator  is  |<A+(0  >,  and  the  intermediate 
state  produced  by  the  Mpller  Out  operator  is  |i/f-(t)  >.  Because  the  Hilbert  Space  can 
be  divided  into  two  orthogonal  spaces:  the  bounded  and  unbounded  spaces,  every  state 
vector  |i/^(t)  >  that  belongs  to  the  unbounded  (scattered)  space  can  be  produced  from  either 
\^reactant  >  Of  product  >  usiug  the  Mpllcr  In  and  Mpller  out  operators  respectively.  Setting 
the  two  equations  grouped  in  Eq.  2.4  equal  to  each  other  and  using  the  property  that  Mpller 
operators  are  isometric,  the  product  state  vector  can  be  solved  in  terms  of  the  reactant  state 
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vector  as  shown  below: 


(2.6) 


where  S  is  the  S  operator  [6,  8,  9,  16].  The  S  operator,  which  is  a  unitary  operator,  contains 
all  of  the  information  regarding  the  scattering  event.  Thus  the  probability  that  a  reactant 
makes  a  transition  to  a  product  is  given  by  Pi^f  =  \  <  4'product\S \^reactant  >  1^.  In  order  to 
determine  the  matrix  elements  of  the  S  operator,  the  channel  packet  method  is  used. 


2.2  Channel  Packet  Method 

In  order  to  determine  the  matrix  elements  of  the  S  operator,  a  basis  must  be  chosen 
in  order  to  represent  the  operator.  The  asymptotic  Hamiltonian  can  be  separated  into 
a  Hamiltonian  that  describes  the  internal  motion  of  the  reactant  or  product  and  into  a 
Hamiltonian  that  describes  the  relative  momentum  of  the  reactants  or  products  with  the 
internal  states.  If  y  is  a  quantum  number  that  describes  all  the  internal  quantum  states  of 
the  asymptotic  Hamiltonian  and  \ky  >  specifies  the  relative  momentum,  then  the  vector 
\ky,  y  >  spans  the  vector  which  is  known  as  the  y  arrangement  channel,  and  this  vector  is 
an  eigenvector  of  the  asymptotic  Hamiltonian.  When  the  Mpller  operators  act  on  this  basis, 
the  \ky,y  >  is  transformed  into  a  vector  labelled  \ky,y+  >.  Because  of  the  intertwining 
relation  of  the  Mpller  Operators  {CiZ+H^  =  hCi\),  it  can  be  shown  that  the  transformed 
basis,  \ky,  y±  >,  is  an  eigenvector  of  the  full  Hamiltonian.  The  reactant  or  product  state 

expanded  as  a  linear  combination  of  the  \ky,  y  >  basis  is 

00 

>=  J'  f?±(ky)lky,y  >  (2.7) 

—  00 

where  rj+iky)  are  the  expansion  coefficients.  When  the  reactant  or  product  state  is 
transformed  into  the  Mpller  state  via  the  Mpller  operators,  the  expansion  coefficients  of 
the  Mpller  state  with  respect  to  the  \ky,  y+  >  are  the  same  as  the  expansion  coefficients  for 
the  reactant  or  product  state  in  Eq.  2.7. 
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In  order  to  find  the  s-matrix  elements  [8,  9,  11,  16,  17],  the  Fourier  transform  of  the 


time  evolution  of  the  reactant  M0ller  state,  labelled  >,  is  taken: 


00 

-f 


IAl(E)>=  I  e^\ifjl> 


(2.8) 


Substituting  the  expansion  coefficients  into  Eq.  2.8,  using  the  fact  that  \ky,y+  >  is  an 
eigenvector  of  the  full  Hamiltonian,  taking  the  projection  of  Eq.  2.8  onto  the  product  Mpller 
state,  and  using  the  orthogonality  relations  of  \ky,y+  >,  Eq.  2.8  reduces  to 

In/u 


<iffl\Al(E)  >= 


QkyWkyl) 


(2.9) 


W_(+k;,)u,(+ky)S^2 

y'’  7 

+r]l(+k',)r]^(-ky)S'(^ 

y  ' 

+ril(-k;,)r^,(+ky)S^:i 

y 

+C(-k;,)77+(-k,)S!';^^_,^] 

where  the  +  and  -  sign  represent  asymptotic  states  with  positive  or  negative  momentum. 
Erom  Eq.  2.9,  we  see  that  careful  consideration  of  the  expansion  coefficients  allows  one 
to  isolate  a  particular  s-matrix  element.  Eor  example,  if  the  reactants  were  to  only  have 
momentum  in  the  direction  of  the  interaction  potential  and  if  the  products  were  to  only 
have  momentum  away  from  the  interaction  potential,  then  only  the  would  survive 

because  r]+(-ky)  =  rj*_(+ky)  =  0.  If  the  expansion  coefficients  were  chosen  with  care,  the 
s-matrix  elements  are 


cr'r 

^±k'  ,,±ky 
y'  ^ 


h\\ky\\ky\)-2 

2nnT]*_(±ky)j]+(±ky) 


<  >k^'\Al(E)  > 


(2.10) 


where  <  if/-  |A+(£')  >  is  the  Eourier  transform  of  <  if/l  >,  and  <  ij/Z.  \e^\ij/\  >  is 

the  correlation  function,  Cyyit).  The  correlation  function  is  the  overlap  between  evolving 
reactant  Mpller  state  at  each  point  in  time  with  a  product  state  that  is  stationary  in  time  for  a 
particular  arrangement  channel.  The  correlation  function  is  a  time  signal  that  depends  upon 
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the  accuracy  of  the  propagation  scheme  of  the  reactant  Mpller  state  [8,  9,  11,  16,  17].  In 
Eq.  2.10,  the  correlation  function  is  used  as  a  tool  to  calculate  s-matrix  elements;  however, 
because  the  correlation  function  is  akin  to  a  signal  in  time  that  contains  information  about 
the  scattering  event,  other  time  signal-processing  tools  can  be  performed  on  the  correlation 
function  to  hopefully  determine  information  about  the  scattering  event.  One  such  signal¬ 
processing  technique  involves  the  Wigner  Distribution  Function  [21]. 
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III.  Wave  Packet  Propagation 


In  order  to  produce  correlation  functions  that  can  be  used  to  generate  Wigner 
Distribution  Functions,  the  propagation  of  the  reactant  Mpller  state  must  be  accurate. 
The  Split  Operator  Equation  will  be  used  to  propagate  the  reactant  Mpller  state,  and  the 
accuracy  of  the  Split  Operator  Method  will  be  verified  with  the  propagation  of  a  wave 
packet  in  free  space  and  a  coherent  state  in  a  quantum  harmonic  oscillator. 

3.1  Split  Operator  Theory 

The  time  evolution  of  a  quantum  system  is  governed  by  Schrdinger’s  Equation: 

|(/r(?)  >=  =  0)  >  .  (3.1) 

One  method  to  finding  an  analytical  solution  to  >  consists  of  constructing  the 
Hamiltonian  and  then  constructing  the  appropriate  transformation  matrix  that  diagonalizes 
the  Hamiltonian.  The  transformation  operator  is  applied  to  Eq.  3.1,  and  >  is 
transformed  to  the  eigenvector  basis  of  the  Hamiltonian.  Eor  the  scattering  problem,  we 
wish  to  instead  express  the  Hamiltonian  in  terms  of  the  kinetic  and  potential  operators, 

I{f+V)t 

W(t)  >=  e  h  =  0)  >  (3.2) 

where  f  is  the  kinetic  energy  operator  and  V  is  the  potential  operator.  The  Taylor  Series 
expansion  of  Eq.  3.2  is  given  by: 

|(A(ao  >=  [1  -  — (T  +  y)  -  — (T^  +  ry  +  yr  +  y)  +  ...]m  =  O)  >  .  (3.3) 

To  third  order,  this  equation  is  equivalent  to  the  Taylor  Series  expansion  of 

,  .  -iVAt  -if  At  -iVAl  .  ,  ^ 

|i/^(At)  >=  e  e  n  e  in  |^(^  =  0)  >  .  (3.4) 

Eq.  3.4  is  called  the  Split  Operator  Equation  [7-9,  11,  17].  The  potential  operator  is 
diagonalized  in  the  position  basis,  and  the  kinetic  operator  is  diagonalized  in  the  momentum 
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basis.  In  order  to  study  the  time  evolution  of  a  wave  packet,  we  need  to  apply  a 
transformation  to  the  momentum  basis  and  back  to  the  position  basis.  The  transformation 
between  position  and  momentum  basis  is  the  inverse  Fourier  transform. 

3.1.1  Fourier  Transform  and  Discrete  Fourier  Transform. 

The  Fourier  transform  belongs  to  a  class  of  transformations  that  transforms  a  vector  in 
one  basis  to  a  vector  in  another  basis.  For  example,  the  inverse  Fourier  transform  transforms 
the  function  ^(x)  =<  x\i//  >  from  the  coordinate  representation  of  the  vector  \tff  >  into 
the  momentum  representation  of  |i/^  >  which  is  =<  k\i//  >.  In  order  to  accomplish 
a  transformation  between  one  representation  of  a  vector  into  another  representation,  the 
projection  of  the  one  basis  onto  the  other  basis  must  be  known.  For  example,  to  obtain  ik(k), 
we  use  completeness  of  the  coordinate  vector  to  show  that  <  k\i//  >=<  k\  f  |x  ><  x\dx\il/  >= 
f  <  k\x  ><  x\iIj  >  dx.  The  projection  of  the  coordinate  vector  onto  the  momentum  vector, 
<  k\x  >,  is  and  from  properties  of  linear  algebra  <  x\k  >=<  k\x  >*=  The  Fourier 
transform,  !F,  is  defined  to  be  a  transformation  from  one  basis  to  another  such  that  the 
projection  of  the  continuous  basis  onto  the  other  continuous  basis  is  and  the  inverse 
Fourier  transform,!F“^  will  have  the  projection  equal  to  In  terms  of  functions,  the 
Fourier  transform  looks  like 

00 

ijjik)  =  J"  i]/ix)e~‘^^dx  (3.5) 

—  00 
00 

ikix)  =  ^J  ^{k)e‘'^"dk.  (3.6) 

— oo 

The  factor  f-  arises  from  the  fact  that  the  Fourier  transform  and  the  inverse  Fourier 
transform  are  not  symmetric.  Without  the  two  pi  factor,  TT~^\}k{xy\  +  iA(x),  and  the 
result  would  be  off  by  a  constant  factor  [3]. 

The  Fourier  transform  assumes  that  the  basis  vectors  are  continuous;  however,  due  to 
the  limits  of  a  computer,  a  discrete  set  of  points,  N,  are  sampled.  As  a  result,  the  integral 
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must  be  approximated  by  a  discrete  sum,  and  the  Fourier  transform  and  its  inverse  are  now: 

N-l 

i//(k)  =  Ax'^ij/(xj)e‘''^^^  (3.7) 

j=o 

^  N-l 

=  ^  (3.8) 

The  discrete  Fourier  transform  requires  on  the  order  of  calculations.  For  large  arrays 
with  lengths  such  as  10^,  this  can  take  up  to  two  weeks  of  CPU  time.  In  order  to  create 
a  more  efficient  algorithm  to  calculate  the  discrete  Fourier  transform,  the  Fast  Fourier 

Transform  was  developed  (FFT).  The  idea  behind  the  FFT  is  that  a  discrete  Fourier 

transform  of  length  N  can  be  rewritten  as  the  sum  of  two  discrete  Fourier  transforms,  each 
of  length  Y-  This  process  can  be  done  recursively  such  that  each  component  is  further 
divided  into  two  discrete  Fourier  transforms  each  with  half  the  length.  If  the  original 
number  of  sampled  points,  N,  is  an  integer  power  of  two,  then  the  discrete  Fourier  transform 
can  be  reduced  to  a  series  of  transforms  of  length  one.  Using  bit  reversal  and  combining 
pairs,  the  discrete  Fourier  transform  can  be  calculated  on  the  order  of  Nlog2N  calculations 
which  reduces  the  CPU  time  to  30  seconds  for  an  array  of  length  10^.  If  the  number  of 
sampled  points  is  not  an  integer  power  of  two,  then  Press  [3]  recommends  that  data  set  is 
at  a  minimum  padded  with  zeros  until  the  length  of  the  data  set  is  a  power  of  two. 

3.1.2  Fast  Fourier  Transform. 

Due  to  the  sampling  theorem  [3],  the  grid  sizes  of  conjugate  variables  are  linked  in  an 
FFT.  The  coordinate  grid  will  consist  of  n  evenly  spaced  points  ranging  from  a  minimum 
value,  Xmin,  to  a  maximum  value  x^ax-  The  spacing  of  the  coordinate  grid.  Ax  is  simply, 

.  ^max  ^niin 

Ax  = - ; — .  (3.9) 

n-l 

The  maximum  and  minimum  values  of  k  which  depend  upon  the  sampling  interval.  Ax,  are 
given  by  the  Nyquist  critical  momentum,  and  the  evenly  spaced  interval  of  the  momentum 
grid  also  depends  upon  the  resolution  of  the  coordinate  grid  and  the  number  of  sampled 
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points.  The  momentum  grid  is  defined  as: 


^maxlmin  —  — 


Ak 


2n 

2n 


nAx 


(3.10) 

(3.11) 


As  shown  in  Eq.  3.10,  increasing  the  resolution  in  the  coordinate  system  (decreasing 
Ax)  decreases  the  resolution  of  the  momentum  grid  and  the  bandwidth  of  momentum 
that  can  be  sampled.  The  increments  between  data  points  in  each  grid  system  must  be 
capable  of  resolving  the  features  of  the  wavefunction  in  the  basis  defined  by  the  grid. 
Similarly,  the  range  of  momentum  for  a  particular  wavefunction  must  be  within  the  bounds 
established  by  the  Nyquist  critical  momentum.  Otherwise,  distortions  will  appear  in  the 
EFT  because  the  momentum  that  fall  outside  the  range  in  Eq.  3.10  will  be  falsely  translated 
into  the  domain.  This  phenomenon  is  called  aliasing  [3,  7].  To  avoid  aliasing,  the  initial 
conditions  of  the  wavefunction  must  have  momentum  that  falls  to  zero  within  machine 
error  near  the  maximum  and  minimum  values  of  momentum,  or  the  wavefunction  must 
contain  at  least  two  points  per  cycle  of  the  highest  momentum  value  present.  Although 
the  Fourier  transform  between  the  momentum  and  coordinate  representations  of  a  function 
were  focused  on,  the  Fourier  transform  can  also  transform  a  time  signal  into  a  frequency 
spectrum,  and  the  same  considerations  applies.  The  EFT  used  in  this  project  was  FFTPackS 
which  is  a  Fortran  subroutine  library  writen  by  Paul  Swarztrauber  and  Richard  Valent.  The 
FFT  package  is  independent  of  the  grid  that  the  Fourier  transform  is  taking  place  on,  so  it 
is  up  to  the  user  to  properly  analyze  the  resulting  grid. 


3.2  Free  Particle  Propagation 

A  particle  in  free  space  experiences  no  potential  by  definition  [4] .  As  a  result,  there 
are  no  boundary  conditions  to  restrict  the  energy  eigenvalues  of  the  Hamiltonian,  and  the 
particle  is  expected  to  have  a  continuous  range  of  energy  values.  Although  the  plane  wave 
is  a  solution  to  Schrdinger’s  equation  to  a  particle  in  free  space,  the  plane  wave  does  not 
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belong  to  the  Hilbert  space  because  the  plane  wave  is  not  a  square  integrable  function. 
This  means  that  the  plane  wave  is  not  a  physically  realizable  quantum  state.  Instead  the 
superposition  of  plane  waves  which  is  square  integrable  is  used  as  a  solution  to  Schrdinger’s 
Equation  [19].  Mathematically,  a  this  superposition  of  plane  waves  which  is  called  a  wave 
packet  has  the  following  general  form: 

00 

<  x\i]/{t)  >=  N  f  c{k)e''^^dk  (3.12) 


where  N  is  the  normalization  factor  for  the  wave  packet,  c{k)  =<  k\il/{t  =  0)  >  are 
the  expansion  coefficients,  and  is  the  plane  wave.  One  set  of  expansion  coefficients, 
c{k)  =  -  kof  -  iKxo]  where  6  is  the  spread  of  Gaussian  wave  packet 

projected  onto  the  coordinate  representation  (the  uncertainty  in  x  or  the  width  of  the 
Gaussian),  xq  is  the  initial  displacement  of  the  Gaussian,  and  k^  is  the  initial  displacement 
of  the  Gaussian  wave  packet,  creates  a  special  type  of  wave  packet  called  the  Gaussian 
wave  packet.  The  evolution  of  a  free  particle  that  was  initially  in  a  Gaussian  state  will  have 
the  following  analytical  solution: 


+  iko{x  -  Xo) 


ik^th 


O  1  iht  1  4(52  '  2li 

tk(x,t)  =  (2n6^)  ?(1  +  ;^)  ^Exp[ - ^ - - - —]. 

2p6  1  +  ^ 


(3.13) 


The  Gaussian  wave  packet  projected  onto  the  momentum  representation  which  is  the 
inverse  Fourier  transform  of  Eq.  3.13  is 


25^  1  o  0  ik^ht 

ikik,  t)  =  ( — YExp[-6^{k  -  ko)  -  ikoXo  -  — — ] 

TT  2/i 


(3.14) 


The  propagation  of  a  Gaussian  wave  packet  is  known  analytically  as  described  by  Eq.  3.13. 
Because  the  particle  is  in  free  space,  the  potential  is  assumed  to  be  zero  for  all  values  of 
position,  and  the  potential  operator  reduces  to  the  zero  operator.  As  a  result,  the  kinetic 
operator  now  commutes  with  the  potential  operator,  and  the  wavefunction  determined  by 
the  split  operator  method,  Eq.  3.4,  reduces  to 


|(/^(dO  >=  e  \4t{t  =  0)  >  . 


(3.15) 
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which  is  exact.  Eq.  3.13  and  Eq.  3.15  will  produce  the  same  result  regardless  of  the  step 
size  in  time.  The  kinetic  operator  is  diagonalized  in  the  momentum  representation.  The 
values  in  a  diagonalized  matrix  are  the  eigenvalues  of  the  associated  operator  that  the  matrix 
represents.  Because  \if/(t  =  0)  >  is  initially  projected  onto  the  coordinate  representation, 
an  EET  from  the  coordinate  representation  of  the  initial  wavefunction  to  the  momentum 
representation  is  required.  Because 

F{A)\(P  >=  F{a)\(p  >  (3.16) 

where  F  is  a  function  of  the  operator  A,  a  is  the  eigenvalue  of  the  operator,  and  \(p  >  is 
the  eigenvector  of  the  operator,  \if/(6t)  >  is  a  diagonalized  matrix.  In  order  to  find  the 
wavefunction  at  a  particular  time  >,  the  process  is  iterated  with  \il/{n6t)  >  becoming 
the  initial  condition  for  |</^((n  +  1)50  >  where  n  is  an  integer.  Generally,  the  coordinate 
representation  is  preferred,  so  another  Eourier  transform  is  required  to  transform  > 
into  the  coordinate  representation. 

The  parameters  of  the  initial  Gaussian  wave  that  was  propagated  along  with  the 
propagation  parameters  and  grid  sizes  are  in  Table  3.L  The  purpose  of  the  free  particle 
propagation  was  to  test  the  propagation  and  EET  codes  for  a  simple  case,  so  the  parameters 
were  chosen  for  simplicity  and  have  no  resemblance  to  any  physical  phenomenon.  The 
wave  packet  was  propagated  for  2  atomic  units  and  compared  to  the  analytical  results. 
Eig.  3.1(a)  shows  the  initial  wave  packet  and  Eig.  3.1(b)  shows  the  analytical  and 
computation  results.  As  evident  in  Eig.  3.1(b)  ,  the  analytical  and  computation  results 
overlap.  Eig.  3.2  shows  the  relative  error  between  the  propagated  wavepacket  and  the 
analytical  wavepacket. 


16 


Table  3.1:  Gaussian  Wave  Paeket  Parameters 


variable 

value 

meaning 

^time 

2048 

number  of  time  grid  points 

Xo 

0 

initial  reaetant  eoordinate  displaeement 

ko 

0 

initial  reaetant  momentum  displaeement 

6 

1 

initial  spread  of  reaetant  in  x 

tmin 

0 

minimum  time 

^max 

2 

maximum  time 

Pceition  (A.U)  Position  (A.U) 


(a)  The  initial  Gaussian  wavepacket 


(b)  The  propagated  wavepacket 


Figure  3.1:  The  norm  of  the  initial  Gaussian  wavepaeket,  with  parameters  deseribed  in 
table  3.1,  propagated  for  a  time  of  2  atomie  units  in  eoordinate  spaee.  The  norm  of  the 
evolved  Gaussian  wavepaeket,  blue  squares,  is  plotted  on  the  same  figure  as  the  norm  of 
the  analytieal  solution,  red. 
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0.12 


Figure  3.2:  The  relative  error  between  the  propagated  wavepaeket  and  the  analytieal 
wavepaeket. 


3.3  Propagation  of  a  Coherent  State  of  a  Quantum  Harmonic  Oscillator 

A  coherent  state  of  the  harmonie  oscillator  is  defined  as  the  quantum  state  where  the 
expectation  value  of  observed  quantities  equals  the  classical  value  of  the  same  observable 
for  all  time  in  the  limit  of  large  quantum  numbers  [4].  For  example,  the  momentum  and 
position  of  a  classical  particle  in  a  harmonic  potential  oscillate  between  a  minimum  and 
maximum  value.  The  expectation  values  of  the  momentum  and  position  operators  of  a 
coherent  state  will  match  these  values  for  all  time.  Using  the  displacement,  raising,  and 
lowering  operators,  it  can  be  shown  that  the  displacement  operator  acting  on  the  ground 
state  of  the  quantum  harmonic  oscillator  is  a  coherent  state.  The  displacement  operator 
simply  boosts  and  translates  the  ground  state,  which  is  a  Gaussian  with  parameters  that 
depend  upon  the  harmonic  potential.  If  <  x  >  is  the  expectation  value  of  the  coordinate  x, 
<  P  >  is  the  expectation  value  of  momentum  P,  co  is  the  frequency  of  the  oscillator,  m  is 
the  mass  of  the  particle,  cr  is  a  complex  number,  and  0  is  a  phase  factor,  then  the  coherent 
state  for  the  harmonic  oscillator  is  defined  as 


= 


mo)^  1 
nn 
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(3.17) 
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Because  the  ground  state  of  the  harmonic  oscillator  is  stationary  state,  the  time  evolution 
of  the  coherent  state  is  straightforward: 

<x>(t)=  ^Re[a(t)] 

<  p  >  (t)  =  ^2mh(x>lm\a{t)\ 


Ax  = 

/  — 

Y  2mci) 

AP  = 

/  mhco 

V  2 

Oait)  = 

4 

We  see  that  when  lal  >>  1,  the  expectation  values  of  the  position  and  momentum  operators 
equal  the  classical  trajectories.  A  particle  in  a  classical  harmonic  trajectory  will  oscillate 
between  a  maximum  position  and  a  minimum  position  where  momentum  is  zero  at  these 
positions.  The  center  of  the  oscillation  will  have  the  largest  momentum.  From  Eq.  3.18, 
the  width  and  amplitude  of  the  wave  packet  does  not  vary  with  time.  Physically,  this  is 
due  to  the  fact  that  spreading  of  the  Gaussian  wave  packet  is  overcome  by  the  squeezing 
of  the  harmonic  potential.  This  squeezing  of  the  wave  packet  from  the  harmonic  potential 
results  from  the  tendency  to  move  away  from  the  edges  of  the  potential  where  the  potential 
is  large  [4].  This  feature  of  the  coherent  state  is  very  valuable  because  the  Gaussian  wave 
packet  can  remain  localized  to  the  harmonic  oscillator  throughout  the  entire  oscillation  and 
because  the  norm  of  the  Gaussian  wave  packet  {ij/  *  if/)  translates  back  and  forth  without 
changing  shape. 

Because  the  propagation  of  a  coherent  state  of  a  quantum  harmonic  oscillator  is  known 
analytically  for  all  time,  the  solution  arrived  at  using  the  split  operator  method  can  be 
compared  with  analytical  results,  Eq.  3.18.  A  coherent  state  solution  is  a  displaced  ground 
state  as  discussed  previously,  and  the  displaced  ground  state,  which  is  a  Gaussian  wave 
packet,  will  oscillate  with  expectation  values  matching  classical  trajectories.  Eor  our  code. 
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a  Gaussian  wave  packet  was  placed  at  the  amplitude  of  oscillation,  and  the  resulting  wave 
packet  was  recorded  at  various  times. 

Because  the  potential  operator  and  the  kinetic  operator  do  not  commute,  a  set  of 
basis  vectors  cannot  simultaneously  diagonalize  both  the  kinetic  and  potential  operators. 
The  basis  for  which  the  potential  operator  is  diagonal  is  the  coordinate  basis.  Because 
the  initial  condition  is  given  in  the  coordinate  representation,  we  can  apply  Eq.  3.16 
directly  to  the  initial  condition.  Because  the  kinetic  operator  is  diagonalized  in  the 
momentum  representation,  an  EFT  is  used  to  transform  e'^\ij/{t  =  0)  >  into  the  momentum 
representation.  After  applying  Eq.  3.16,  another  EFT  is  used  to  transform  the  result 
back  into  the  coordinate  representation  since  the  potential  operator  is  diagonalized  in  the 
coordinate  representation.  The  evolution  of  the  initial  wavefunction  at  a  small  time  step 
later,  |i/^(At)  >,  is  known  after  applying  Eq.  3.16  one  last  time.  The  code  to  propagate  the 
initial  wave  packet  now  looks  like 

if/(x,At)  =  e—r[e—r~\e—i//(x,0)]]  (3.19) 

In  order  to  determine  the  wavefunction  at  a  particular  time  t,  Eq.  3.19  is  calculated 
iteratively  [7-9].  Because  the  amplitude  and  width  of  the  coherent  state  remains  constant 
and  because  the  coherent  state  can  oscillate  on  a  grid  that  fits  within  the  computational 
window,  the  evolution  of  the  coherent  state  can  be  observed  for  large  times.  Small  errors 
that  occur  in  a  single  iteration  of  the  split  operator  method  will  grow  and  after  many 
iterations  can  become  large.  We  can  use  this  concept  to  determine  the  if  the  size  of  the 
time  step  used  in  the  split  operator  was  sufficient  for  the  split  operator  to  approximate.  It 
was  shown  by  Eacy  [7]  and  verified  here  that  a  time  step  of  approximately  At  =  0.01  atomic 
units  is  sufficient  for  the  split  operator  approximation. 

The  parameters  of  the  initial  coherent  state,  the  propagation  parameters,  the  grid  sizes, 
and  the  potential  parameters  are  shown  in  Table  3.2.  The  purpose  of  the  coherent  state 
propagation  was  to  test  the  propagation  code,  so  the  parameters  were  again  chosen  for 
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simplicity  and  have  no  resemblance  to  any  physical  system.  Because  the  coherent  state 
follows  the  classical  trajectory  of  a  particle  in  the  same  potential,  the  initial  coherent  state 
was  set  up  at  the  maximum  amplitude  of  oscillation.  The  coherent  state  was  propagated  for 
t  =  50.25:;r  atomic  units  and  compared  to  the  analytical  results  for  the  corresponding  times. 
The  frequency  of  the  oscillator  is  1  atomic  unit.  Fig.  3.3(a)  shows  the  initial  coherent  state 
and  Fig.  3.3(b)  shows  the  analytical  and  computation  results.  As  evident  in  Fig.  3.3(b) ,  the 
analytical  and  computation  results  overlap. 


Table  3.2:  Coherent  State  Parameters 


variable 

value 

meaning 

^time 

4096 

number  of  time  grid  points 

^min 

0 

minimum  time 

^max 

50.25;7r  maximum  time 

<  Xo> 

5 

initial  coordinate  displacement 

<  ko  > 

0 

initial  momentum  displacement 

6 

0.707 

initial  spread  of  reactant  in  x 

CO 

1 

frequency  of  oscillator 

1 

mass  of  particle 
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Figure  3.3:  The  Coherent  State  of  quantum  harmonie  oseillator.  The  potential  parameters 
and  the  initial  eondition  of  the  eoherent  state  are  found  in  table  3.2.  The  initial  eoherent 
state  whieh  is  a  displaeed  Gaussian  wavepaeket  is  shown  in  the  figure  on  the  left.  The 
eoherent  state  is  propagated  to  a  time  of  50.25;m>  atomie  units  in  eoordinate  spaee  in  the 
figure  on  the  right.  The  numerieally  eomputed  eoherent  state,  blue  squares,  is  plotted  in  the 
same  figure  on  the  right  as  the  analytieal  solution,  red.  Both  figures  overlap. 
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Figure  3.4:  The  relative  error  between  the  propagated  eoherent  state  at  time  t  =  50.257rm 
atomie  units  and  the  analytieal  eoherent  state. 
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IV.  Square  Well 


4.1  Square  Well 

The  square  well  potential  allows  for  both  bounded  states  {E  <  0)  and  seattering  states 
(E  >  0)  [4,  5];  however,  we  will  only  eonsider  seattering  states.  The  square  well  potential 
has  three  regions  on  the  position  grid.  In  the  first  and  third  regions,  the  potential  is  assumed 
to  be  zero;  however,  in  the  seeond  region,  the  potential  is  a  eonstant  value,  V,  over  a  length 
of  L.  Eaeh  of  these  regions  assumes  a  solution  to  the  wavefunetion  in  the  form  of 

il/(x)  =  (4.1) 

where  i  labels  the  region  and  A  and  B  are  eonstants  in  the  region.  Eq  4.1  is  the 
superposition  of  an  ineoming  plane  wave  with  an  outgoing  plane  wave  for  eaeh  value  of 
momentum.  An  ineoming  wave  is  defined  as  a  plane  wave  heading  towards  the  square  well 
and  the  outgoing  wave  is  defined  as  a  plane  heading  away  from  the  square  potential. 

At  the  boundaries  of  the  square  well,  the  potential  jumps  from  the  eonstant  potential 
to  the  zero  potential;  however,  the  wavefunetion  must  be  eontinuous  and  the  first  derivative 
of  the  wavefunetion  must  also  be  eontinuous.  Applying  the  eontinuity  eonditions  at  both 
boundaries  of  the  well  boundaries,  we  have  a  system  of  four  equations  with  six  unknowns. 
Solving  the  eoeffieients  of  the  first  region  in  terms  of  the  eoeffieients  in  the  seeond  region, 
and  solving  the  seeond  region’s  eoeffieients  in  terms  of  the  third  regions  eoeffieients,  the 
first  region  ean  be  solved  in  terms  of  the  third  region  for  a  given  energy.  The  eonneetion 
between  the  first  and  third  region  ean  be  aeeomplished  with  a  2  x  2  transformation  matrix 
for  eaeh  value  of  momentum,  k.  The  first  and  third  regions  are  akin  to  the  asymptotie 
regions  deseribed  in  the  seattering  theory  seetion,  and  the  seeond  region  is  equivalent  to 
the  interaetion  region.  Instead  of  a  matrix  that  transforms  the  third  region  into  the  first 
region,  we  seek  a  matrix,  the  s-matrix,  that  eonneets  the  ineoming  waves  to  the  outgoing 
waves.  The  outgoing  waves  of  the  first  and  third  region  ean  be  solved  using  algebra  in 
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terms  of  the  outgoing  waves  of  the  first  and  third  region  for  a  given  energy.  If  we  assume 
that  there  are  no  incoming  waves  in  the  third  region,  then  the  outgoing  waves,  Bi  and  A3  , 
can  be  expressed  in  terms  the  incoming  wave,  Ai  using  only  two  s-matrix  elements.  The 
s-matrix  element,  that  relates  A3  with  Ai  is  related  to  the  transmittance  of  Ai  through 
the  square  well  while  the  s-matrix  element,  Ski,ki,  that  relates  Bi  with  Ai  is  related  to  the 
reflectance  of  Ai  through  the  square  well.  The  transmission  coefficient  is 


T(E)  =  \Sk,M\ 


1 


1  +  ^2m{E  +  V)) 


(4.2) 


4E{E+V)  h 

where  L  is  the  length  of  the  well,  m  is  the  mass  of  the  particle,  E  is  energy,  and  V  is 
the  magnitude  of  the  well  depth.  Eq.  4.2  provides  an  analytical  solution  to  compare  the 
computationally  produced  transmission  coefficient.  Because  the  computationally  produced 
transmission  coefficient  involved  the  channel  packet  method,  Eq.  4.2  was  used  to  verify  the 
accuracy  of  the  correlation  function  produced. 

The  sinusoidal  term  in  Eq.  4.2  creates  a  series  of  oscillating  peaks  and  valleys  between 
the  values  of  one  and  zero.  The  location  of  the  peaks  in  the  transmission  coefficient  align 
with  the  eigenvalues  of  the  infinite  square  well  potential.  The  eigenvalues  of  the  infinite 
square  well  potential  are  En  =  where  n  is  the  eigenvalue,  so  the  location  of  the 
peaks  is  given  by  E  =  E„-V.  The  location  of  the  valleys  of  the  transmission  coefficient  can 
be  approximated  as  the  average  between  two  consecutive  peaks.  The  difference  between 
the  peak  and  valley  describes  the  resonance  of  the  square  well.  The  larger  the  difference 
between  the  peak  and  valley,  the  greater  the  resonance.  The  degree  of  resonance  in  the 
peaks  of  the  transmission  coefficient  is  quantized  by  a  parameter  value,  /3  which  depends 
upon  the  characteristics  of  the  square  well  as  shown  below  [24]: 

LV2my 


h 


(4.3) 


Re-expressing  Eq.  4.2  in  terms  of  jS  and  the  eigenvalues  of  the  infinite  square  well,  we  have 

1 


T{E)  =  \Sk,M\  = 


1  + 


4nrfi(En-V) 


sin^(^  V2m(E  -I-  V)) 


(4.4) 
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We  see  that  /3  influences  the  amplitude  of  the  sinusoidal  term.  When  the  argument  of 
the  sine  function  in  Eq.  4.4  is  zero,  the  transmission  coefficient  is  a  maximum,  but  when 
the  argument  of  the  sine  function  is  |  radians,  the  transmission  coefficient  is  a  minimum. 
The  amplitude  of  this  minimum  is  determined  by  the  amplitude  of  the  sine  function  which 
is  •  As  yS  increases,  the  difference  between  the  maximum  and  minimum  values 

increases.  From  Eq.  4.3,  we  see  that  increasing  the  depth  of  the  well  increases  the 
(5  parameter  or  resonance  of  the  transmission  coefficient  which  was  produced  from  the 
correlation  function. 

In  order  to  calculate  the  correlation  function,  a  wave  packet  was  propagated  through  a 
one-dimensional  finite  square  well  potential  in  order  to  calculate  s-matrix  elements  which 
determine  the  transmission  coefficient.  To  summarize  the  propagation  scheme: 

1.  Set  up  the  initial  wavefunction,  |i/r(0)  >,  in  coordinate  space:  <  .r|^(0)  > 

2.  Operate  on  the  wavefunction  with  e~^ 

3.  Transform  the  resulting  wavefunction  to  momentum  space  via  an  inverse  EFT 

4.  Operate  on  the  wavefunction  with  e'^ 

5.  Transform  the  resulting  wavefunction  to  coordinate  space  via  a  EFT 

6.  Operate  on  the  resulting  wavefunction  with  e~^ 

7.  Repeat  items  2-6  for  each  At  step  until  <  x\il/(t)  >  is  obtained 

The  resolution  of  the  time  grid  must  be  taken  into  consideration  in  order  for  the  above 
propagation  scheme  to  work.  The  split-operator  theory  is  only  accurate  with  the  Taylor 
Series  expansion  of  Schrdinger’s  equation  up  to  third  order  assuming  that  At  is  small.  As 
shown  with  the  harmonic  oscillator,  a  time  step  of  At  «  0.01  is  sufficient  for  propagation. 
Similarly,  the  resolution  of  the  x  grid  is  important  not  only  because  of  Eq.  3.10  but  because 
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of  the  boundary  of  the  square  potential.  For  a  square  potential,  the  potential  function  can 
be  thought  of  as  a  step  function  at  the  edges  of  the  potential.  If  Is.x  was  too  large,  then 
the  distance  between  two  points  on  the  x  grid  would  be  large,  and  the  potential  would  no 
longer  drop  off  like  a  step  function.  Instead,  the  square  well  or  square  barrier  would  become 
a  trapezoidal  potential  which  has  different  properties  than  the  square  potential,  and  the  code 
would  not  match  the  analytical  results. 

4.2  Absorbing  Boundary  Conditions 

Another  consideration  that  needs  to  be  taken  into  account  are  the  boundaries  of  the 
computational  window  [7,  11].  In  Fortran,  if  a  function  extends  passed  the  window,  it 
loops  back  to  the  other  side.  Consequently,  a  propagating  wavefunction  with  positive 
momentum  that  reaches  the  maximum  x  grid  coordinate  will  loop  back  to  the  minimum 
X  grid  coordinate  and  add  to  the  component  of  the  wavefunction  that  was  present  at  Xmin- 
This  issue  has  not  been  a  problem  with  previous  potentials  because  we  were  not  concerned 
with  either  propagating  the  wavefunction  for  long  periods  or  the  wavefunction  was  trapped 
oscillating  within  the  window  for  infinity;  however,  because  calculating  s-matrix  elements 
and  Wigner  Distribution  Functions  requires  the  correlation  function  to  die  off,  the  wave 
packet  must  be  propagated  for  long  periods  of  time.  One  fix  to  the  issue  would  be  to 
extend  the  grid  to  larger  values  of  x.  Unfortunately,  this  solution  comes  with  the  trade-off 
in  computation  time.  In  order  to  be  able  to  keep  the  resolution  of  the  wavefunction,  the 
number  of  sampled  points  must  increase  as  the  size  of  the  grid  increases. 

An  alternative  solution  would  be  to  impose  absorbing  boundary  conditions.  Absorbing 
boundary  conditions  force  the  wavefunction  at  the  boundaries  to  decay  to  zero  in  order  to 
minimize  artificial  contributions  to  the  correlation  function  when  the  wavefunction  loops 
to  the  other  side  [7,  11].  To  accomplish  this,  an  imaginary  potential  term  localized  at 
the  edges  of  the  computation  window  will  be  added  to  the  interaction  potential  such  that 
y  =  yint  +  iVabc  where  Vabc  is  the  absorbing  boundary  potential  and  V  is  the  interaction 
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potential.  When  the  potential  operator  is  exponentiated  in  Eq.  3.4  the  imaginary  component 
becomes  real,  and  if  the  potential  has  eigenvalues  that  are  negative  when  x  is  positive  or 
positive  when  x  is  negative,  the  wavefunction  will  be  dampened.  Although  many  potentials 
are  available  as  candidates,  we  do  not  want  to  drastically  dampen  the  wavefunction  on  one 
extreme  or  slowly  dampen  the  wavefunction  on  the  other  extreme  such  that  the  correlation 
function  is  altered.  As  a  result,  we  choose  the  following  form 

Vabc(x)  =  Vabse~''  ^  +  Vabse^'-  ^  (4.5) 

where  Vabs  is  the  strength  of  the  absorbing  boundary  potential  and  is  the  spread  of  the 
absorbing  boundary  condition.  Eq.  4.5  is  the  sum  of  two  Gaussian  functions  centered  at  the 
edges  of  computational  grid. 

4.3  Reactant  and  Product  States 

The  transmission  coefficient  for  a  finite  square  has  analytical  solution  according 
Eq.  4.2  which  was  compared  to  numerical  results.  In  order  to  determine  the  transmission 
coefficient,  the  correlation  function  between  a  reactant  Mpller  state  and  a  product  Mpller 
state  for  the  finite  square  well  must  be  calculated.  The  reactant  state  was  a  Gaussian  wave 
packet  with  positive  momentum  placed  to  the  left  of  the  square  well.  The  product  state  was 
a  Gaussian  wave  packet  with  positive  momentum  placed  to  the  right  of  the  square  well. 
As  mentioned  previously,  the  channel  packet  method  relies  on  intelligently  choosing  the 
expansion  coefficients  of  the  momentum  representation  of  the  reactant  and  product  states  in 
order  to  isolate  and  solve  for  a  particular  s-matrix  element.  Because  we  are  interested  in  the 
transmission  coefficient,  we  are  interested  in  the  wavefunction  that  propagated  through  the 
square  potential.  Eor  a  reactant  state  placed  to  the  left  of  the  well  and  a  product  state  placed 
to  the  right  of  the  well,  both  states  must  have  non-zero  expansion  coefficients  for  only 
positive  values  of  momentum.  Because  the  reactant  and  product  states  are  Gaussian  wave 
packets,  the  expansion  coefficients,  which  is  a  Gaussian  wave  packet  in  the  momentum 
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representation,  are  known  analytieally  from  Eq.  3.14.  Beeause  the  reaetant  and  produet 
states  were  plaeed  outside  the  well  where  the  potential  is  zero,  the  reaetant  Mpller  and 
produet  Mpller  states  are  the  same  as  the  initial  reaetant  and  produet  state  respeetively.  A 
Gaussian  wave  paeket  propagated  to  infinity  under  zero  potential  and  then  baek  to  its  initial 
loeation  under  zero  potential  will  be  the  same  as  the  initial  Gaussian  as  verified  by  Eq.  3.13. 
Consequently,  \il/ reactant i product  >=  l<A±  >>  and  no  numerieal  eomputations  had  to  be  made  to 
produee  the  Mpller  states. 

The  parameters  of  the  reaetant  state  and  produet  states,  the  propagation  parameters, 
and  grid  sizes  are  in  Table  4.1.  The  parameters  of  the  reaetant  and  produet  states  were 
intended  to  mateh  the  parameters  of  the  H2  moleeule;  however,  the  propagation  times 
neeessary  to  produee  a  eorrelation  funetion  that  eaptured  the  full  seattering  event  required 
millions  of  data  points.  Beeause  the  purpose  of  the  square  well  was  to  provide  eonfidenee 
in  the  aeeuraey  of  the  WDEs  and  to  gain  insight  into  the  affeet  of  potential  depths  (or 
heights)  on  the  features  of  the  WDE,  some  of  the  H2  parameters  were  redueed  in  order  to 
reduee  the  number  of  time  data  points  to  a  manageable  number  to  save  on  eomputation 
time.  A  lighter  partiele  with  large  energies  takes  less  time  to  travel  a  given  distanee  than  a 
heavier  partiele  with  lower  energies.  If  the  resolution  of  the  time  grid  has  to  remain  fixed  at 
approximately  At  =  0.01,  then  a  shorter  range  in  time  required  less  grid  points  to  maintain 
the  resolution.  As  a  result,  a  fourth  of  the  redueed  mass  of  H2  was  used  instead  of  918.076 
atomie  mass  units,  and  the  average  energy  of  the  moleeule  was  3.486  Hartree  energy  whieh 
eorresponds  to  94.844  eV.  The  well  depths  of  the  square  well  were  based  upon  the  well 
depth  of  the  ground  state  of  H2  whieh  is  0.166  Hartree.  The  potential  was  then  varied 
to  determine  the  impaet  of  large  resonant  struetures  in  the  transmission  ooefReient  on  the 
WDE.  The  minimum  well  depth  was  0.166  Hartree  while  the  maximum  well  depth  was 
16.611  Hartree  (100  times  larger  than  the  minimum  well  depth).  Eig.  4.1  shows  the  initial 
set  up  for  the  reaetant  and  produet  states  with  a  well  depth  of  1L628  Hartree.  Eig.  4.3(a) 
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shows  the  momentum  representation  of  the  reactant  state  which  is  the  same  as  the  product 
state.  As  shown  in  Fig.  4.3(b)  ,  is  very  small  dXk  <  which  is  necessary  in  order  to 
solve  for  in  Eq.  2.9. 

y'  ’ 


Table  4.1:  Square  Well  and  Square  Barrier  Parameters 


variable 

square  well  value 

square  barrier  value 

meaning 

^max 

95 

450 

maximum  coordinate 

^min 

-95 

-450 

minimum  coordinate 

^max 

180 

1250 

maximum  time 

^min 

0 

0 

minimum  time 

Xq 

-2.8 

-2.8 

initial  reactant  coordinate  displacement 

k{) 

40 

40 

initial  reactant  momentum  displacement 

6 

0.075 

0.075 

initial  spread  of  reactant  in  x 

xpo 

2.8 

2.8 

product  state  coordinate  displacement 

kpo 

40 

40 

product  state  momentum  displacement 

6p 

0.075 

0.075 

initial  spread  of  product  state  in  x 

P 

229.519 

229.519 

reduced  mass 

a 

4 

4 

well  length 

^abs 

1 

1 

strength  of  absorbing  boundary 

Vd 

2 

2 

spread  of  absorbing  boundary 

n 

8192 

32768 

number  of  spatial  grid  points 

^time 

8192 

131072 

number  of  time  grid  points 

Hr 

8192 

32768 

number  of  tau  grid  points 

ds 

1 

4 

number  of  points  down  sampled 
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Figure  4.1:  The  initial  setup  to  the  square  well  problem.  The  square  well  of  width  4  A.U.  is 
eentered  about  x  =  0.  For  this  example,  the  depth  of  the  square  well  is  1 1 .628  Hartree.  The 
reaetant  state  is  the  blue  Gaussian  wavepaeket  to  the  left  of  the  square  well  and  the  produet 
state  is  the  red  Gaussian  wavepaeket  to  the  right  at  the  square  well.  The  parameters  of  the 
reaetant  and  produet  states  are  found  in  table  4.2 
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Figure  4.2:  The  absorbing  boundary  eonditions  are  shown  as  the  green  Gaussians  at  the 
edge  of  the  eoordinate  grid  for  the  same  well  shown  in  Fig.  4. 1 .  The  parameters  of  the 
reaetant  and  produet  states  are  found  in  table  4.2.  The  absorbing  boundary  eondition  is 
shown  for  the  positive  eoordinate  values;  however,  there  is  a  symmetrie  Gaussian  at  the 
edge  of  the  grid  for  the  negative  eoordinate  values. 
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(a)  The  momentum  representation  of  the  reactant  (b)  A  zoomed  in  graphic  of  the  momentum 
state  representation  of  the  reactant  state 

Figure  4.3:  The  momentum  representation  of  the  reaetant  state  assoeiated  with  the 
parameters  found  in  Table  4.2.  The  expansion  eoeffieients  are  all  positive  valued  and  at 
k  =  0,  the  expansion  eoeffieients  are  eomputationally  equivalent  to  zero 

As  shown  in  Table  4. 1 ,  the  propagation  times  for  the  square  well  problem  was  different 
than  the  propagation  times  for  the  square  barrier.  The  reason  for  this  is  that  a  partiele 
in  the  square  well  gains  kinetie  energy  while  the  partiele  in  square  barrier  looses  kinetie 
energy.  As  a  result,  the  partiele  takes  a  longer  time  to  propagate  through  the  potential  of 
a  square  barrier  than  a  square  well.  If  the  resolution  of  the  time  grid  has  to  remain  fixed 
at  approximately  At  =  0.01,  more  grid  points  are  needed  to  speeify  the  time  grid  for  the 
square  barrier  eompared  to  the  square  well  potential.  Similarly,  the  range  of  the  eoordinate 
grid  inereased  for  the  square  barrier  potential  as  eompared  to  the  square  well  potential  sinee 
a  wave  paeket  with  a  longer  propagation  time  travels  more  distanee,  so  the  number  of  grid 
points  for  the  eoordinate  grid  also  inereased.  An  example  of  the  partiele  propagated  for 
180  atomie  units  through  a  square  well  is  shown  in  Fig.  4.4.  We  ean  see  in  Fig.  4.4  that  part 
of  the  original  wave  paeket  was  refleeted,  part  of  the  wave  paeket  was  transmitted  through, 
and  although  very  small  eompared  to  the  rest  of  the  wavefunetion,  part  of  the  wavefunetion 
exists  in  the  well  region. 
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Figure  4.4:  The  propagated  reactant  state  shown  in  Fig.  4.1  with  a  square  well  depth  of 
8.305F'//.  Part  of  the  wavefunction  transmitted  through  the  well,  part  of  the  wavefunction 
reflected  away  from  the  product  state,  and  part  of  the  wavefunction  still  oscillates  inside 
the  well. 


4.4  Correlation  Function 

The  correlation  function,  Cyyit)  =<  >,  is  the  projection  of  the  evolving 

reactant  Mpller  state  onto  the  product  Mpller  state.  Because  the  Mpller  states  are  the  same 
as  the  initial  states,  the  correlation  function  for  the  square  potential  is 


oo 


(4.6) 


— oo 


where  il/product{x,to)  is  the  product  Mpller  state  in  the  coordinate  representation  at 
a  particular  time  to  and  il/reactant(x,t)  is  the  reactant  Mpller  state  in  the  coordinate 
representation  as  a  function  of  time.  Eq.  4.6  shows  the  computational  approach  to 
calculating  correlation  functions  for  the  square  potential.  The  product  Mpller  state 
remained  fixed  as  the  reactant  Mpller  state  was  then  propagated  through  the  square  well. 
The  correlation  function  or  overlap  between  the  propagating  reactant  Mpller  state  with  the 
fixed  product  Mpller  state  was  computed  for  each  time  step  iteration  of  the  split  operator 
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method.  At  each  time  step,  the  product  of  the  reactant  Mpller  state  with  the  conjugate  of 
the  product  Mpller  state  was  taken  at  each  point  on  the  coordinate  grid.  The  products  were 
summed  over  the  coordinate  grid  in  order  to  determine  a  scalar  value  at  that  particulate 
time.  These  scalar  values  produce  a  signal  known  as  the  correlation  function. 

Because  the  signal  depends  upon  the  overlap  between  the  product  and  reactant  Mpller 
states,  the  signal  is  expected  to  be  very  close  to  zero  as  the  initial  states  are  localized 
on  either  side  of  the  square  potential  (the  initial  wavefunctions  die  off  to  the  same  order  of 
magnitude  as  machine  error  at  the  edges  of  the  potential).  As  the  reactant  state  is  propagated 
through  the  square  well  or  barrier,  the  signal  increases  and  peaks  as  the  wavefunction  passes 
through  the  product  state.  Due  to  reflections  at  the  boundaries  of  the  square  potential,  some 
of  the  amplitude  of  the  wave  packet  will  oscillate  in  the  region  of  the  potential  before  being 
transmitted  through  or  reflected  backwards  out  of  the  potential.  Consequently,  there  is 
asymmetry  in  the  correlation  function  as  components  of  the  reactant  Mpller  state  overlap 
with  the  product  Mpller  state  at  later  times.  The  wavefunction  must  be  propagated  for  a 
long  time  in  order  to  ensure  that  oscillating  components  of  the  wavefunction  inside  the 
potential  are  captured  and  the  full  scattering  event  is  collected  in  the  signal. 

As  the  reactant  state  was  propagated  through  the  potential,  the  correlation  function 
was  calculate  at  each  position.  Fig.  4.5(a)  and  Fig.  4.5(b)  shows  an  example  of  the  real 
and  imaginary  parts  of  a  correlation  function  for  a  potential  well  depth  of  11.628  Hartree. 
In  Fig.  4.5,  the  correlation  function  is  zero  for  approximately  the  first  20  atomic  units  in 
time.  This  makes  sense  since  the  reactant  and  product  Mpller  states  do  not  overlap  each 
other.  At  approximately  20  atomic  units  in  time,  the  higher  energy  plane  waves  that  made 
up  the  Gaussian  wave  packet  travelled  through  the  square  well  potential  and  entered  the 
region  where  the  product  Mpller  state  was  placed.  As  time  increases,  a  greater  number  of 
plane  waves  with  a  larger  spread  in  energy  enter  the  region  where  the  product  Mpller  state 
is  localized  increasing  the  amplitude  of  overlap.  As  the  transmitted  wave  packet  continues 
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to  propagate,  the  higher  energy  plane  wave  components  leave  region  where  the  product 
Mpller  state  is  placed,  and  the  norm  of  the  correlation  function  decreases.  At  approximately 
30  atomic  units  of  time,  the  correlation  function  is  approximately  zero,  as  most  of  the 
transmitted  wave  packet  passed  through  the  localized  region  of  the  product  Mpller  state. 
However,  during  this  time,  a  percentage  of  the  amplitude  of  the  wave  packet  was  reflected 
at  the  boundary  of  the  square  well.  This  reflected  wave  packet  travels  the  square  well  and 
either  reflects  or  transmits  each  time  the  wave  packet  reaches  the  boundary  of  the  square 
well.  After  the  smaller  in  amplitude  wave  packet  traversed  the  length  of  the  square  well 
three  times  (is  at  the  positive  value  boundary  a  second  time),  a  percentage  of  the  wave 
packet  amplitude  transmits  through  the  square  well  and  into  the  product  Mpller  state  region 
for  a  second  time.  This  occurs  at  approximately  40  atomic  units.  Because  a  percentage  of 
the  wave  packet  either  transmits  or  reflects  at  each  boundary,  the  amplitude  of  the  wave 
packet  that  overlaps  with  the  product  Mpller  state  for  a  second  time  is  significantly  smaller 
compared  to  the  initial  fraction  of  the  wave  packet  amplitude  that  transmitted  through  the 
well.  However,  the  higher  energy  plane  wave  components  that  make  up  the  wave  packet 
will  reach  the  product  Mpller  state  before  the  lower  energy  plane  wave  components  do, 
and  the  correlation  function  will  have  another  peak.  This  process  continues  to  repeat,  and 
in  Fig.  4.5  a  third  peak  is  still  noticeable  before  the  amplitude  of  wave  packet  is  too  small 
compared  to  the  scale  of  the  initial  transmitted  wave  packet.  The  correlation  function  is 
approximately  zero  for  times  greater  than  80  atomic  units. 
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Figure  4.5:  Correlation  Function  Plots  from  the  overlap  between  the  reactant  state  and 
product  state  with  parameters  given  in  Table  4.2  and  a  potential  of  1 1 .628  Hartree 
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Figure  4.6:  The  norm  of  the  correlation  function  shown  in  Fig.  4.5. 


4.5  Down  Sampling 

According  to  the  sampling  theorem,  at  least  three  points  are  needed  to  be  sampled 
for  each  peak  or  valley  in  order  to  capture  a  particular  feature  of  a  continuous  function 
[3].  Because  the  split  operator  method  requires  the  time  grid  to  be  sampled  with  a  very 
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small  interval  between  each  point,  At,  the  correlation  function  may  contain  a  very  large 
number  of  points.  Although  a  greater  number  of  sampled  points  more  accurately  reflects 
the  continuous  function,  there  is  a  trade  off  with  computational  time.  In  some  cases 
where  a  large  number  of  points  described  each  feature  of  the  correlation  function,  it  may 
be  necessary  to  down  sample  the  correlation  function  in  order  to  significantly  reduce  the 
computational  time  of  the  correlation  function.  Down  sampling  is  the  elimination  of  evenly 
spaced  data  points  from  the  total  set  of  data  points.  In  this  paper,  a  down  sample  of  one 
means  that  every  point  was  used  while  a  down  sample  of  two  signifies  that  every  other  point 
was  used. 

4.6  Transmission  Coefficient 

As  discussed  previously,  the  correlation  function  can  be  used  to  determine  the  s-matrix 
elements.  In  order  to  calculate  s-matrix  elements,  the  FFT  of  the  correlation  function 
was  taken  with  respect  to  the  time  variable  and  normalized  by  the  expansion  coefficients 
and  momentum  as  given  in  Eq.  2.10.  Energy  can  be  related  back  to  the  momentum  of 
a  free  particle  in  order  to  arrive  at  the  s-matrix  elements  as  a  function  of  energy  rather 
than  momentum.  Eor  these  calculations,  only  \S +k+k\'^  was  calculated.  Other  s-matrix 
elements  such  as  |5  -k+k\^  can  be  computed  if  the  initial  set-up  was  different.  Eor  example, 
if  the  product  Mpller  state  was  instead  located  to  the  left  of  the  square  well  with  negative 
momentum,  IS-k+k?'  can  be  isolated  in  Eq.  2.9.  Einally,  the  transmission  coefficient, r(£'), 
is  IS+k+kiE)]^.  Comparing  the  numerically  calculated  transmission  coefficient  with  the 
analytical  solution  provided  confidence  that  the  computed  correlation  function  was  correct. 
It  is  important  to  note  that  when  the  Eourier  transform  of  the  correlation  function,  c(co),  was 
divided  by  the  expansion  coefficients  in  the  region  where  both  the  expansion  coefficients 
and  c(m)  were  smaller  than  the  machine  error,  the  s-matrix  elements  blew  up.  As  a  result, 
there  is  a  range  of  energy  values  that  correspond  to  a  realistic  transmission  coefficient  (i.e 
the  range  of  values  for  which  the  T{E)  is  less  than  one). 
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With  the  correlation  functions  calculated,  the  transmission  coefficients  and  scattering 
WDFs  were  calculated.  The  square  well  did  not  require  any  down  sampling;  however,  the 
square  barrier  did.  According  to  Table  4.1,  every  fourth  point  of  the  correlation  function 
was  used.  Fig.  4.6  shows  the  norm  squared  of  the  Fourier  transform  of  the  correlation 
function  used  in  Fig.  4.5  with  the  expansion  coefficients  plotted  on  the  same  scale,  and 
Fig.  4.6  is  an  example  of  the  transmission  coefficient  calculated  from  the  correlation 
function  in  Fig.  4.5.  As  shown  in  Fig.  4.6,  the  transmission  coefficient  is  only  valid  for 
a  range  of  energy  values.  Outside  the  range  of  values,  the  transmission  coefficient  blows 
up  because  the  values  of  c(oj)  and  of  the  expansion  coefficients  are  smaller  than  the  error 
in  producing  the  correlation  function.  Note  that  because  we  are  in  atomic  units,  energy  and 
frequency  are  interchangeable  since  E  =  hco.  Also  note  that  the  transmission  coefficient 
oscillates  between  one  and  a  value  less  than  one.  When  T{E)  =  1,  the  well  is  completely 
transparent  to  plane  waves  of  energy,  E.  The  /3  parameter  describes  how  resonant  each  of 
these  peaks  are.  The  greater  the  p  parameter,  the  sharper  and  more  spread  out  these  peaks 
in  the  transmission  coefficient  become. 
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Figure  4.7 :  [The  norm  squared  of  the  Fourier  transform  of  the  correlation  function  is  plotted 
as  a  function  of  frequency.  The  correlation  function  was  the  overlap  between  the  reactant 
state  and  product  state  with  parameters  given  in  Table  4.2  and  a  potential  of  1 1 .628  Hartree. 
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Figure  4.8:  The  transmission  coefficient  as  a  function  of  energy.  The  computationally 
produced  transmission  coefficient,  blue  squares,  is  compared  to  the  analytically  determined 
transmission  coefficient,  red  fine.  The  error  is  due  to  the  resolution  in  the  coordinate 
grid.  Because  the  coordinate  grid  is  sampled  at  discrete  points,  the  finite  square  well  is 
not  perfectly  square.  There  is  a  slight  angle  to  the  edges  of  the  well  making  it  slightly 
trapezoidal. 


4.7  Autocorrelation  Function 

The  correlation  function  was  then  used  to  produce  the  Wigner  Distribution  Function. 
The  scattering  WDF  is  the  Fourier  transform  of  the  non-stationary  autocorrelation  function 
of  the  correlation  function  which  is  expressed  in  terms  of  the  time,  t,  and  a  time  lag,  r. 
As  shown  in  Eq.  4.7,  the  autocorrelation  function  is  c*{t  +  \)c{t  -  \).  The  autocorrelation 
function  is  the  product  of  the  correlation  function  shifted  in  the  positive  time  direction 
with  its  own  complex  conjugate  shifted  by  the  same  amount  except  in  the  negative  time 
direction.  Theoretically,  this  product  is  summed  for  every  shift  ^  ranging  from  -oo  to 
oo;  however,  computationally  this  can’t  be  done.  Instead,  the  shifts  can  only  be  done  in 
discrete  points  over  a  finite  range.  The  range  of  r  must  be  large  enough  such  that  the 
shifted  correlation  functions  have  very  little  overlap  at  the  maximum  values  of  r.  At  the 
maximum  or  minimum  values  of  r,  the  product  c*{t  +  ^)c(t  -  ^)  needs  to  be  nearly  zero  for 
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all  values  of  t.  At  r  =  t^ax  -  tmin,  the  distance  between  c*(t  +  and  c(t  -  is  the  length  of 
the  signal,  and  the  product  will  be  c*(t  +  ^)c(t  -  for  all  values  of  time.  This  implies  that 
the  appropriate  range  for  r  is  -(t^ax  -  tmin)  +  2At  <  r  <  (t,„ax  -  t^m)- 

Another  limitation  is  that  the  correlation  function  is  an  array  of  length  n,,  the  number 
of  time  points.  Any  shift  has  to  be  done  such  that  c(t,  -  =  c(tj).  In  other  words,  the  shifted 

correlation  function  needs  to  still  align  with  the  same  sampled  points  on  the  time  grid.  This 
implies  that  At  =  2{tmax  -  tmin)-  In  addition,  the  correlation  array  is  of  finite  length  with 
each  element  of  the  array  corresponding  to  a  particular  point  on  the  time  grid.  If  an  indexed 
value  is  shifted  to  a  time  outside  the  time  grid,  then  the  autocorrelation  function  is  assumed 
to  be  zero  at  that  point.  In  order  to  determine  if  the  integrand  was  properly  calculated, 
the  autocorrelation  function  at  a  time  lag  of  r  =  0  should  be  equivalent  to  the  magnitude 
squared  of  the  correlation  function,  |c(t)|^. 

The  real  and  imaginary  parts  of  the  autocorrelation  function  associated  with  the 
correlation  function  in  Fig.  4.5  are  plotted  as  contours  in  Fig.  4.9.  The  real  part  can  be 
seen  to  exhibit  an  even  function  centered  about  t  =  0  atomic  units  while  the  imaginary  part 
is  an  odd  function  centered  about  r  =  0.  Thus,  the  scattering  WDF  must  be  a  real  function. 
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(a)  The  Real  Part  of  c*(t  +  ^)c(t  -  ^)  (b)  The  Real  Part  of  c*(t  +  §)c(f  -  f ) 

Figure  4.9:  Autocorrelation  Function  Plot  from  the  overlap  between  the  reactant  state  and 
product  state  with  parameters  given  in  Table  4.1  and  a  potential  of  11.628  Hartree.  The 
real  part  of  the  autocorrelation  function  is  an  even  function  centered  about  r  =  0  while  the 
imaginary  part  of  the  autocorrelation  function  is  an  odd  function  centered  about  t  =  0. 
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4.8  Phase  Shift 


In  order  to  produce  the  scattering  WDF,  the  FFT  of  the  autocorrelation  function  was 
taken  with  respect  to  the  time  delay.  Because  the  autocorrelation  function  has  an  even 
real  part  and  an  odd  imaginary  part,  symmetry  analysis  tells  us  that  the  scattering  WDF 
must  be  a  purely  real  function  which  is  a  property  that  is  always  true  for  the  scattering 
WDF.  Unfortunately,  if  FFT  of  the  autocorrelation  function  was  taken  without  further 
consideration,  the  scattering  WDF  would  have  an  imaginary  component.  The  reason 
for  this  is  that  the  FFT  package  assumes  that  the  minimum  value  of  t  is  zero  instead  of 
Tmin  =  -(tmax  “  tmin)  +  2At.  As  a  rcsult,  the  cvcnncss  and  oddness  of  a  function  disappears. 
In  order  to  properly  account  for  the  shift  in  the  time  grid,  a  phase  shift  must  be  added  to 
the  result  of  the  FFT.  The  phase  shift  ,which  can  be  deduced  by  applying  a  time  shift  to 
the  inverse  Fourier  transform  (Eq.  3.5),  is  and  when  the  phase  shift  is  taken  into 

account,  the  imaginary  component  of  the  WDF  is  approximately  zero.  The  phase  shift  has 
not  been  an  issue  previously  because  we  have  been  interested  in  the  norm  square  of  any 
FFT  result  where  the  phase  information  is  lost  or  because  the  FFT  result  was  transformed 
back  into  the  original  coordinate  system  using  an  inverse  FFT  which  undoes  the  phase  shift. 

4.9  Wigner  Distribution  Functions 

The  Wigner  Distribution  Function  is  a  distribution  created  by  Eugene  Wigner  in  1932 
that  attempts  to  relate  quantum  mechanics  with  the  classical  phase-space  picture;  however, 
the  Wigner  Distribution  Eunction  (WDE)  can  also  be  found  in  signal-processing  and  optics 
[2,  7,  14,  15,  21].  The  Wigner  Distribution  Eunction  as  found  in  quantum  mechanics  is 
the  Eourier  transform  of  the  density  operator  with  respect  to  a  coordinate  variable.  The 
goal  of  the  WDE  was  to  find  quantum  corrections  to  classical  statistical  mechanics  and 
to  find  a  relation  between  the  probability  densities  in  quantum  with  the  trajectories  in 
phase  space  in  classical  physics.  Due  to  the  uncertainty  principle,  it  is  impossible  to  find 
a  probability  distribution  in  phase  space  such  that  the  probability  distribution  is  positive 
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everywhere  and  that  the  integral  of  the  product  of  the  probability  distribution,  P{x,  y),  and 
of  some  observable,  A(x,p),  that  is  a  function  of  the  position  and  momentum  operators 
produces  the  expectation  value  of  that  operator,  f  f  P(x,  p)A(x,  p)dxdp  =<  A  >.  Thus, 
the  Wigner  Distribution  Function  is  not  a  probability  distribution  function  since  it  can  be 
negative,  but  the  WDF  is  a  quasi-probability  distribution  function.  Although  there  are 
other  quasi-probability  distribution  functions  that  have  been  proposed,  the  application  of 
the  other  quasi-distribution  functions  has  been  found  to  be  limited. 

The  quantum  perspective  of  the  scattering  Wigner  Distribution  Function  is  similar  to 
the  time-signal  point  of  view  [21].  Instead  of  density  operator,  the  integrand  on  the  time 
signal,  or  in  this  case  the  correlation  function,  is  a  non- stationary  autocorrelation  function. 
The  autocorrelation  of  a  signal  is  the  cross-correlation  of  a  signal  with  itself  as  a  function  of 
time  lag.  If  the  signal  is  only  composed  of  one  time  series,  then  the  WDF  can  be  simplified 
to 


Wit,  co) 


+  7,)CYyit  -  “^dT. 


(4.7) 


where  r  is  a  time  lag.  The  scattering  WDF  in  Eq.  4.7  is  the  Fourier  transform  of  the  product 
of  the  correlation  function  shifted  by  an  offset  in  time  with  its  complex  conjugate  shifted 
in  the  opposite  direction  in  time.  Thus,  the  scattering  WDF  in  Eq.  4.7  is  a  quasi-spectral 
distribution  function  of  the  correlation  function;  however,  because  the  correlation  function 
is  the  overlap  between  two  wavefunctions,  the  scattering  WDE  in  Eq.  4.7  can  be  viewed 
as  a  quasi  distribution  of  probabilities  associated  with  a  particular  energy  value.  In  other 
words,  the  scattering  WDE  should  describe  the  flow  of  energy  probabilities  of  a  scattering 
event  as  a  function  of  time.  As  a  result,  a  negative  value  in  the  scattering  WDE  can  be 
interpreted  as  a  decrease  in  probability. 

Some  unique  properties  of  the  scattering  WDE  that  are  of  interest  for  this  work  are 


f  W(t,  a))da)  =  |cyy(t)|^ 
f  W(t,co)dt  =  \cyy(a))\'^ 
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where  Cyyia))  is  the  Fourier  transform  of  the  correlation  function  [2,  14,  15].  Eq.  4.8 
states  that  the  projection  of  the  scattering  WDF  onto  the  time  domain  reproduces  the  norm 
squared  of  the  correlation  function.  Similarly,  the  projection  of  the  scattering  WDF  onto  the 
frequency  domain  reproduces  the  norm  squared  of  the  Fourier  transform  of  the  correlation 
function.  From  the  channel  packet  method,  if  the  Fourier  transform  of  the  correlation  is 
divided  by  the  expansion  coefficients  of  the  reactant  Mpller  state  along  with  some  other 
prefactors,  then  the  s-matrix  elements  are  obtained  as  seen  by  Eq.  2.10.  As  a  result,  the 
scattering  WDF  contains  all  of  the  time  and  energy  information  (frequency  is  related  to 
energy)  regarding  the  scattering  event. 

4.10  WDF  for  Square  Potentials 

The  scattering  WDF  is  the  FFT  of  the  autocorrelation  function  for  every  value  of  time. 
For  a  particular  value  of  t,  the  autocorrelation  function  was  determined  as  a  function  of  r 
which  was  then  sent  into  the  FFT  package.  The  result  was  multiplied  by  a  phase  factor 
and  some  constants  determined  by  Eq.  4.7.  This  process  was  repeated  for  each  point  on 
the  time  grid.  In  order  the  determine  the  validity  of  the  scattering  WDF,  the  integrated 
projections  of  the  scattering  WDF  onto  the  time  and  frequency  domains  should  match  the 
magnitude  squared  fo  the  correlation  function,  |c(0p  and  the  magnitude  squared  of  the 
Fourier  transform  of  the  correlation  function,  |c(a»)p,  respectively  as  shown  in  Eq.  4.8. 
The  integrals  were  simply  sums  multiplied  by  the  respective  grid  interval.  Scattering 
WDFs  were  calculated  for  a  particle  scattering  from  1 1  different  square  well  potentials. 
Each  square  well  potential  had  the  same  length  and  the  same  reactant  and  product  states; 
however,  the  depth  of  the  square  well  varied.  The  purpose  of  producing  these  scattering 
WDFs  for  varying  well  depths  was  to  gain  insight  into  scattering  with  potentials  that  had 
largely  resonant  peaks  in  the  transmission  coefficient  such  as  the  HOCO  scattering  problem. 
Scattering  WDFs  were  also  produced  for  the  finite  square  barrier  with  varying  barrier 
heights.  The  only  changes  required  to  the  program  were  the  sign  on  the  potential  and 
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the  number  of  time  grid  points  used.  The  various  potential  heights  and  associated  (3  values 
are  shown  in  Table  4.2. 


Table  4.2:  Square  Well  and  Square  Barrier  Potentials  Tested 


Well  Depth 

Barrier  Height 

P 

0.16610929 

0.16610929 

l3o  =  34.929 

1.66109294 

1.66109294 

VlO;0o 

3.32218588 

3.32218588 

V^/3o 

4.98327882 

4.98327882 

V30/3o 

6.64437176 

6.64437176 

V4O;0o 

8.30546470 

8.30546470 

V^/3o 

9.96655764 

9.96655764 

11.62765058 

11.62765058 

13.28874352 

13.28874352 

14.94983646 

14.94983646 

V%ySo 

16.61092940 

16.61092940 

10>So 

In  Fig.  4.10,  the  scattering  WDF  is  plotted  as  a  function  of  time  and  energy.  The 
plot  is  a  color  image  of  the  scattering  WDF  associated  with  the  correlation  function  in 
Fig.  4.5.  Although  the  absolute  magnitude  of  the  scattering  WDF  is  not  easily  read  from 
the  color  image  in  Fig.  4. 10,  the  location,  shape,  and  relative  heights  of  certain  features  of 
the  scattering  WDF  are  much  easier  to  see.  In  Fig.  4.10,  we  can  see  a  large,  negatively 
chirped  peak  followed  by  bands  of  oscillating  peaks  and  valleys  at  later  times.  This  chirped 
peak  is  positioned  such  that  at  high  energies  and  short  times,  the  scattering  WDF  has  non¬ 
zero  values  and  at  low  energies  and  long  times,  the  scattering  WDF  also  has  non-trivial 
values.  This  makes  sense  when  one  considers  the  various  components  of  the  propagating 
wave  packet.  The  wave  packet,  which  is  a  superposition  of  plane  waves,  has  the  high  energy 
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plane  waves  traveling  faster  than  the  low  energy  plane  wave  (henee  why  the  wave  paeket 
spreads  in  free  spaee).  The  plane  waves  whieh  get  a  boost  from  the  square  well  will  make  it 
to  the  loeation  of  the  stationary  produet  Mpller  state  mueh  sooner  than  the  slower  traveling 
plane  waves.  As  expeeted,  the  high  energy  eomponents  of  the  seattering  WDF  have  peaked 
values  at  larger  times  as  eompared  to  the  lower  energy  eomponents  of  the  seattering  WDF. 
In  Fig  4.10,  one  ean  also  observe  a  seeondary  band  of  peaks  later  in  time  eorresponding  to 
the  same  energy  values  as  the  large  ehirped  peak.  This  smaller  band  of  non-zero  seattering 
WDF  values  has  a  spread  that  ean  be  deseribed  the  by  same  time-energy  relation  as  the 
large  initial  ehirped  peak  diseussed  later. 
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Figure  4. 10:  A  eolor  image  of  the  Wigner  Distribution  Funetion  for  the  eorrelation  funetion 
ealeulated  with  a  square  well  potential  of  1 1 .628  Hartree  and  reaetant  and  produet  states 
with  parameters  speeified  in  Table  4.2. 


The  projeetions  of  the  seattering  WDF  onto  the  time  and  frequeney  domains 
reproduees  the  norm  squared  of  the  eorrelation  funetion  and  the  norm  squared  of  the  Fourier 
transform  of  the  eorrelation  funetion.  Fig.  4.11  and  Fig.  4.12  shows  the  overlap  between  the 
projeetions  and  the  assoeiated  funetions,  \c{t)\^  and  |c(m)|^.  In  the  figure  with  the  projeetion 
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of  the  WDF  onto  the  frequeney  domain,  there  is  slight  disagreement  between  the  projeetion 
of  the  seattering  WDF  onto  the  frequeney  domain  and  \c{(jof\.  The  seattering  WDF  is  the 
Fourier  transform  of  the  time  delay  present  in  the  autoeorrelation  funetion.  The  interval 
between  sampled  points  of  the  r  grid  is  twiee  the  size  as  the  time  grid, At  =  2At.  Beeause 
the  energy  grid  is  inversely  related  to  the  time  grid,  the  spaeing  between  energy  points  for 
the  projeetion  of  the  seattering  WDF  onto  the  frequeney  domain  (Fourier  transform  of  r)  is 
half  the  spaeing  of  \c{oSf\  (Fourier  transform  of  t).  This  implies  that  there  is  an  additional 
point  for  the  projeetion  of  the  WDF  between  two  points  for  |c(m)^|.  When  the  same  points 
are  eompared  (i.e.  every  other  point  of  the  projeetion  of  the  seattering  WDF),  the  projeetion 
of  the  seattering  WDF  and  |c(m)^|  are  in  agreement. 


Figure  4.11:  The  plot  of  the  projeetion  of  the  seattering  WDF  (blue)  onto  the  time  domain 
for  the  square  well  potential  with  a  depth  of  1 1 .628  Hartree  is  shown.  This  projection  is 
compared  to  |c(0p|  (red) 


With  confidence  that  the  WDFs  are  correct,  we  began  generating  scattering  WDFs 
for  many  square  well  and  square  barrier  cases.  Fig  4.13-Fig  4.20  highlight  features  of  the 
scattering  WDF  that  change  as  the  magnitude  of  the  square  potential  increases.  We  see  the 
same  general  structure  in  these  figures  that  we  saw  in  Fig.4.10.  Each  scattering  WDF  had  a 
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Figure  4.12:  The  plots  of  the  projection  of  the  WDF  (blue)  onto  the  time  and  frequency 
domain  for  the  square  well  depth  of  11.628  Hartree  are  shown.  These  projections  are 
compared  to  |c(0p|  and  |c(m)^|  respectively  (red).  The  slight  disagreement  in  amplitude 
between  the  projection  onto  the  frequency  domain  and  \c{cof'\  results  from  differences  in 
sampled  energy  values 


large  chirped  peak  with  the  same  general  shape  as  Fig.  4.10.  The  difference  between  each 
scattering  WDF  is  that  as  the  depth  of  the  square  well  increased,  the  large  peak  shifted 
to  smaller  values  of  time,  and  the  concavity  of  its  chirp  decreased.  As  the  well  depth 
increases,  the  well  adds  more  kinetic  energy  and  hence  more  momentum,  k  =  jq 

the  evolving  reactant  Mpller  state  as  it  propagates  in  the  interaction  region.  As  a  result,  the 
wavefunction  gets  transmitted  through  the  region  much  faster  and  the  wavefunction  that  is 
first  transmitted  through  the  well  should  arrive  sooner  than  for  a  shallower  well.  In  addition, 
as  the  well  depth  increases,  bands  with  oscillating  peaks  and  valleys  appear  like  in  Fig  4.10. 
The  location  of  these  bands  correspond  to  the  same  moments  in  time  as  the  peaks  in  the 
norm  of  the  correlation  function.  As  the  square  well  depth  increased,  the  number  of  visible 
bands  grew,  the  peak  amplitude  decreased,  and  the  difference  between  the  maximum  and 
minimum  values  of  the  scattering  WDF  decreased. 
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Figure  4.13:  Scattering  WDF  for  square  well  depth  of  0.166  Hartree. 
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Figure  4.14:  Scattering  WDF  for  square  well  depth  of  4.9835  Hartree. 


The  scattering  WDFs  for  the  square  barrier  potentials  have  a  similar  structure  as  the 
scattering  WDFs  for  the  square  well.  One  of  the  differences  is  that  as  the  barrier  height 
increases,  the  momentum  of  the  plane  waves  that  make  up  the  wave  packet  decrease.  By 
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Figure  4.15:  Scattering  WDF  for  square  well  depth  of  13.289  Hartree. 
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Figure  4.16:  Scattering  WDF  for  square  well  depth  of  16.61 1  Hartree. 

the  time  the  barrier  height  reached  13.289,  the  transmitted  wavefunction  was  so  small  that 
the  correlation  function  was  on  the  same  order  of  magnitude  as  the  error.  Thus,  scattering 
WDFs  were  only  calculated  up  to  11.628  Hartree.  One  of  the  major  differences  is  that 
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Figure  4.17:  Scattering  WDF  for  square  barrier  height  of  0.166  Hartree.  The  scale  is  based 
off  of  the  maximum  and  minimum  values  of  the  scattering  WDF. 
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Figure  4.18:  Scattering  WDF  for  square  barrier  height  of  3.322  Hartree.  Note  that  the  scale 
is  the  same  as  Fig.  4.17.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  barrier  height  increases. 


the  bands  of  oscillating  peaks  and  valleys  extend  very  far  out  in  time  as  the  barrier  height 
increases  rather  than  coalescing  together.  Another  major  difference  is  that  as  the  barrier 
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Figure  4.19:  Scattering  WDF  for  square  barrier  height  of  4.983  Hartree.  Note  that  the  scale 
is  the  same  as  Fig.  4.17.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of 
the  scattering  WDF  as  the  barrier  height  increases.  On  this  scale  the  scattering  WDF  is  not 
visible. 


Figure  4.20:  Scattering  WDF  for  square  barrier  height  of  4.983  Hartree  scaled  to  the 
maximum  and  minimum  values  of  the  scattering  WDF.  Although  the  amplitude  of  the 
scattering  WDF  if  Fig.  4.19  cannot  be  seen  with  the  scale  used  in  Fig.  4.17,  the  scattering 
WDF  still  exists  as  shown  in  this  figure  . 
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height  increased,  the  curvature  of  the  bands  grew  compared  to  the  flattening  out  of  the  of 
the  bands  in  the  square  well  as  viewed  along  the  energy  axis,  so  the  width  of  the  scattering 
WDF  along  the  energy  axis  decreased  while  also  shifting  to  higher  energy  levels. 

4.10.1  Interpretation  of  Peaks  and  Valleys. 

Since  the  Wigner  Distribution  Function  can  be  negative,  there  hasn’t  been  a  clear 
interpretation  of  the  WDF  in  the  past.  We  conjecture  that  the  positive  values  of  the 
WDF  correspond  to  presence  of  probability  at  the  location  of  the  product  Mpller  state 
with  a  specific  energy  at  a  certain  time,  the  negative  values  correspond  to  the  depletion  of 
probability  at  the  same  product  Mpller  state  with  a  specific  energy  at  a  certain  time,  and 
the  zero  value  corresponds  to  the  presence  of  probability  not  at  the  product  Mpller  state. 
We  arrive  at  this  conjecture  from  looking  at  the  properties  of  the  WDF  defined  by  Eq.  4.8. 
The  projection  of  the  WDF  onto  the  frequency  domain  is  the  norm  Fourier  transform  of  the 
correlation  function  which  is  related  to  the  transmission  coefficient  by  a  normalization  term 
as  shown  in  Eq.  2.10.  According  to  Eq.  4.8,  each  point  on  |c(a»)p  is  the  sum  of  all  values  of 
the  WDE  at  that  particular  frequency  (energy)  over  all  time.  Positive  values  will  contribute 
to  the  sum  at  that  particular  energy  while  negative  values  will  take  away  from  that  sum. 
If  all  the  valleys  in  the  WDE  are  on  the  same  order  of  magnitude  as  the  peaks  for  a  given 
energy  value,  then  the  sum  will  be  close  to  zero  since  each  valley  will  approximately  cancel 
with  each  peak  in  the  sum  (the  amplitude  of  the  peaks  and  valleys  for  each  band  diminishes 
over  time).  This  sum  will  correspond  to  a  near  zero  value  for  |c(a»)p  and  hence  a  very  low 
transmission  coefficient  at  that  associated  frequency(energy). 

As  shown  in  Eig.  4.21  (a  square  well  with  a  well  depth  of  4.984  Hartree),  the 
amplitudes  of  the  valleys  are  not  on  the  same  order  of  magnitude  as  the  amplitudes  of 
the  peaks.  Eor  a  given  energy  value,  |c(m)|^  will  be  large,  and  although  there  are  peaks  and 
valleys  in  the  transmission  coefficient,  the  valleys  in  the  transmission  coefficient  are  near 
one  as  shown  in  Eig.  4.23.  In  contrast  as  shown  in  Eig.  4.22  (a  square  well  with  a  well 
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depth  of  16.611  Hartree),  the  amplitudes  of  the  valleys  are  much  closer  to  the  amplitudes 
of  the  peaks.  As  such,  there  are  energy  values  such  as  3  Hartree  where  T{E  =  3)  is  going 
to  be  much  smaller  compared  to  the  same  transmission  coefficient  for  a  square  well  depth 
of  4.984  Hartree.  The  square  well  with  a  well  depth  of  16.611  Hartree  is  said  to  have 
a  much  deeper  resonance  (much  larger  p  value)  than  the  square  well  of  depth  of  4.984 
Hartree.  Because  the  positive  values  of  the  WDF  contribute  to  the  s-matrix  element  at  a 
given  energy,  we  say  that  the  positive  values  of  the  WDF  must  correspond  to  the  presence 
of  probability  for  a  given  value  of  energy  at  a  given  value  of  time.  On  the  other  hand,  the 
negative  values  of  the  WDF  must  correspond  to  the  loss  of  probability  for  a  given  value  of 
energy  at  a  particular  value  of  time.  Because  the  total  probability  must  be  conserved,  this 
depletion  of  probability  must  show  up  elsewhere  in  the  system  for  a  given  time,  at  some 
other  moment  in  time  if  the  same  probability  shows  up  later  or  earlier  in  time  at  the  same 
location  (product  Mpller  state),  or  both.  The  question  of  where  this  probability  went  will 
be  examined  next. 
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Figure  4.21:  WDF  for  square  well  depth  of  4.984  Hartree.  The  axis  is  scaled  to  the 
maximum  and  minimum  values  of  the  WDF. 
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Figure  4.22:  WDF  for  square  well  depth  of  16.611  Hartree.  The  axis  is  scaled  to  the 
maximum  and  minimum  values  of  the  WDF 
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Figure  4.23:  WDF  for  square  well  depth  of  4.984  Hartree. 
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Figure  4.24:  WDF  for  square  well  depth  of  4.984  Hartree. 


4.10.1.1  Location  of  Lost  Probability. 

If  the  probability  amplitude  for  a  given  energy  is  not  transmitted  at  a  given  time,  then 
it  must  either  be  transmitted  at  a  later  time  or  appear  in  the  refieetion  eoeffieient  whieh 
is  IS-k+kl"^.  To  see  if  the  loss  of  probability  ended  up  in  the  refieetion  eoeffieient,  the 
Wigner  Distribution  Funetion  associated  with  the  reflection  coefficient  had  to  be  produced. 
This  was  accomplished  by  placing  the  same  product  Mpller  state  used  in  the  transmission 
coefficient  to  the  left  side  of  the  well  with  the  following  changes  made  to  Table  4.1: 
xpo  =  -2.8  and  kpo  =  40. 

Fig.  4.26  shows  the  reflection  coefficient  and  Fig.  4.27  shows  that  the  sum  of  the 
transmission  coefficient  and  the  reflection  coefficient  is  one  which  gives  us  confidence  that 
the  correlation  function  for  the  new  set  up  is  correct.  The  WDF  associated  with  \S-k+k\^  is 
plotted  for  a  square  well  with  a  depth  of  11.628  Hartree  as  shown  in  Fig.  4.25.  This  WDF 
is  very  similar  in  structure  to  the  transmission  coefficient.  There  is  a  large  central  peak  and 
several  bands  comprised  of  oscillating  peaks  and  valleys.  The  reflection  WDF  has  a  large 
peak  before  the  transmission  coefficient  since  there  will  be  reflection  at  the  left  edge  of  the 
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square  well.  The  remaining  bands  in  the  reflection  WDF  align  in  time  with  the  bands  in 
the  transmission  coefficient.  This  makes  sense  for  a  deep  square  well  because  as  the  time 
it  takes  the  wave  packet  to  traverse  the  square  well  (t  =  ^2(^v)  ^  the  magnitude 

of  the  well  depth  and  L  is  the  length  of  the  well)  is  very  small  compared  to  the  time  it  takes 
for  the  wave  packet  to  traverse  the  zero  potential  region  (t  =  where  Xo  is  the  distance 
from  the  well  to  the  center  of  the  peak  value  of  the  initial  states). 

Fig.  4.28  shows  the  amplitude  of  the  Wigner  Distribution  Function  at  a  particular 
energy  as  a  function  of  time  (energy  slice)  for  both  the  transmission  and  the  reflection 
WDF.  The  peak  of  the  transmission  coefficient  may  correspond  to  a  valley  in  the  reflection 
coefficient;  however,  this  is  not  true  for  all  peaks  and  valleys.  Fig.  4.29  shows  the  amplitude 
of  the  WDF  at  a  particular  moment  in  time  as  a  function  of  energy  (time  slice)  for  both  the 
transmission  and  reflection  WDF.  The  peaks  and  valleys  of  the  WDF  do  not  align.  In 
fact,  the  frequency  of  peaks  in  the  reflection  WDF  is  double  the  frequency  of  peaks  as 
the  transmission  coefficient.  As  shown  in  Fig.  4.30,  when  the  second  transmission  band 
is  compared  to  the  second  reflection  band  (instead  for  example  as  the  first  transmission 
band  with  the  second  reflection  band  in  Fig.  4.29),  the  peaks  of  the  transmission  WDF 
correspond  to  the  valleys  of  the  reflection  WDF  and  vice  versa.  At  high  energies,  the  peaks 
and  valleys  no  longer  match  up,  but  this  is  due  to  the  chirp.  It  may  be  possible  that  after  the 
initial  reflection  and  transmission,  the  loss  of  probability  in  the  transmission  band  may 
be  found  in  the  reflection  band;  however,  further  study  must  be  done  in  order  to  confirm 
this. 
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Time  (A.U) 

Figure  4.25:  WDF  for  square  well  depth  of  1 1.628  Hartree  and  the  product  state  placed  to 
the  left  of  the  well.  The  correlation  function  and  WDF  are  then  related  to  and  the 

reflection  coefficient. 


Figure  4.26:  Reflection  coefficient  for  square  well  depth  of  1 1.628  Hartree. 
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Figure  4.27:  The  sum  of  the  refleetion  and  transmission  eoeffieient  for  square  well  depth  of 
11.628  Hartree.  The  sum  of  the  refleetion  and  transmission  eoeffleients  is  one  as  expeeted. 


Figure  4.28:  The  amplitude  of  the  transmission  and  refleetion  seattering  WDFs  at  an  energy 
value  of  3.033  Hartree  as  a  funetion  of  time  for  the  square  well  potential  with  a  depth  of 
11.628  Hartree. 
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Figure  4.29:  The  amplitude  of  the  transmission  and  refleetion  seattering  WDFs  at  a  time  of 
32  A.U.  as  a  funetion  of  energy  for  the  square  well  potential  with  a  depth  of  1 1 .628  Hartree. 


Figure  4.30:  The  amplitude  of  the  transmission  and  refleetion  seattering  WDFs  at  a  time 
of  32  A.U.  and  21  A.U.  for  the  square  well  potential  with  a  depth  of  11.628  Hartree. 
The  purpose  of  this  figure  is  to  eompare  the  probability  exehange  between  the  seeond 
transmission  and  refleetion. 


58 


4.10.2  Frequency  of  Arrival. 

Can  the  Wigner  Distribution  Function  reveal  to  us  more  information  about  the 
scattering  event  than  the  correlation  function  or  the  s-matrix  elements?  In  order  to  answer 
this  question,  we  examined  the  Fourier  transform  of  the  amplitude  of  the  WDF  for  a  given 
value  of  energy  over  all  of  time  as  shown  in  Fig.  4.31.  The  slice  in  energy  produces  a 
signal,  and  as  seen  in  Fig.  4.31,  the  Fourier  transform  of  these  signals  shows  the  spectrum 
of  frequencies  present  in  the  signals.  For  the  energy  slice,  the  Fourier  transform  of  this 
signal  will  describe  the  frequency  at  which  probability  was  transmitted  through  the  square 
well.  Hence,  we  will  refer  to  the  Fourier  transform  as  the  frequency  of  arrival  for  a  given 
energy  value.  The  peaks  present  in  Fig.  4.31  is  the  frequency  at  which  a  plane  wave  with 
a  given  energy  oscillates  in  the  well.  This  frequency  is  one  over  the  time  it  takes  for  the 
wave  packet  to  propagate  one  round  trip  in  the  square  well.  The  frequency  is  given  as 


ln^Jl{E  +  V) 

2L^ 


(4.9) 


where  V  is  the  magnitude  of  the  depth  of  the  well.  The  first  peak  corresponds  to  the 
fundamental  frequency  co  =  0.276  A.U.  while  other  peaks  correspond  to  the  odd  harmonics 
of  o).  The  frequency  of  arrival  is  information  that  is  not  found  in  the  correlation  function  of 
the  Fourier  transform  of  the  correlation  function.  As  a  result,  the  WDF  does  provide  more 
information  into  the  scattering  event  than  the  correlation  function  and  s-matrix  elements. 


4.10.3  Negative  Chirp. 

The  shape  of  the  negative  chirp  of  the  WDF  can  be  related  to  the  classical  motion  of 
a  particle.  If  the  distance  from  the  peak  of  the  reactant/product  Mpller  state  to  the  edge  of 
the  square  well  is  Xo,  the  length  of  the  square  well  is  L,  the  velocity  of  the  classical  particle 
in  the  zero  potential  regions  is  v,  and  the  velocity  of  the  classical  particle  in  the  square  well 
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Frequency 

Figure  4.31:  The  figure  is  the  Fourier  transform  of  the  amplitude  of  the  WDF  for  a  square 
well  potential  with  a  depth  of  11.628  Hartree.  The  amplitude  of  the  WDF  was  taken  at 
3.033  Hartree  over  the  entire  time  domain.  The  first  peak  eorresponds  to  the  frequeney  at 
whieh  the  wave  paeket  oseillates  in  the  well  (m),  the  seeond  peak  eorresponds  to  (3m),  the 
third  peak  eorresponds  to  (5m)  and  the  trend  eontinues. 


is  v',  then  the  distanees  the  partiele  travelled  is  given  by 


Xo  =  vAt 
L  =  VAt' 


(4.10) 

(4.11) 


With  E  =  and  E  +  V  =  \iJ.V^,  the  total  time,  t,  ean  be  solved  for  as  a  funetion  of 
energy. 


nr 

t  =  + 


V2(£  +  V) 

where  V  is  the  magnitude  of  the  depth  of  the  square  well. 


) 


(4.12) 
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Figure  4.32:  The  energy-time  relation  of  the  WDF  for  a  square  well  potential  with  a  depth 
of  1 1.628  Hartree.  Eq.  4.12  is  plotted  on  top  of  the  scattering  WDF  as  the  black  line. 
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V.  H  +  H2  and  OH  +  CO  Scattering 


After  producing  scattering  WDF  for  the  finite  square  well  and  square  barrier 
potentials,  scattering  WDFs  were  produced  for  the  H  +  H2  reaction.  Three  scattering  WDFs 
were  produced  for  the  correlation  functions  between  a  reactant  Mpller  state  in  the  ground 
vibrational  state  and  a  product  Mpller  state  in  each  of  the  first  three  internal  vibrational 
energy  states  of  H  +  H2  collision.  The  correlation  function  for  this  scattering  reaction  was 
previously  calculated  by  Weeks  and  Calfas  [11],  and  code  was  developed  to  import  the 
correlation  functions.  In  order  to  calculate  the  correlation  function.  Weeks  and  Calfas  [11] 
used  the  channel  packet  method.  The  reactant  state  was  centered  on  the  Liu-Siegbahn- 
Truhlar-Horowitz  (LSTH)  Potential  Energy  Surface  (PES)  and  coupling  was  assumed  to 
occur  among  the  first  three  surfaces  [12,  13].  Unlike  the  square  potential,  the  reactant 
Mpller  and  product  states  will  not  be  the  same  as  the  initial  reactant  and  product  states 
while  on  these  surfaces.  The  initial  reactant  and  product  states  were  chosen  such  that  the 
computed  Mpller  states  were  localized  in  the  interaction  potential.  In  order  to  compute 
the  Mpller  states,  the  initial  reactant  and  product  states  were  propagated  to  ’’infinity” 
under  the  asymptotic  Hamiltonian  and  then  back  to  zero  under  the  full  Hamiltonian.  The 
correlation  function  is  computed  by  propagating  the  reactant  Mpller  state  forward  in  time 
from  t  =  0  to  tmax  and  then  backwards  in  time  from  t  =  0  to  with  the  aide  of 
absorbing  boundary  conditions.  Because  scattering  can  occur  into  any  of  the  first  three 
internal  vibrational  energy  states,  three  product  Mpller  states  had  to  be  computed,  and  the 
correlation  function  was  computed  for  each  of  these  product  Mpller  states  as  the  reactant 
Mpller  state  evolved  in  time.  As  a  result,  three  scattering  WDEs  were  computed  for  each 
of  the  correlation  functions  in  the  same  way  that  the  scattering  WDEs  were  calculated 
for  the  square  potentials.  The  purpose  of  creating  scattering  WDEs  for  the  H  +  H2  was 
to  verify  the  accuracy  of  the  code  to  import  correlation  functions  and  produce  scattering 


62 


WDFs  from  these  imported  functions  by  comparing  the  scattering  WDF  of  the  ground 
internal  vibrational  arrangement  channel  with  the  scattering  WDF  independently  calculated 
for  the  same  arrangement  channel.  Another  purpose  was  to  produce  scattering  WDFs  for 
the  coupling  between  the  ground  state  with  the  first  two  excited  internal  vibrational  states 
which  was  never  accomplished  before.  The  correlation  function  for  a  single  arrangement 


Figure  5.1:  The  OH{v  =  0)  +  CD  ^  OCO{v  =  0)  +  //  PES. 

channel  of  the  OH{v  =  0)  +  CD  OCO{y  =  0)  +  H  channel  collision  [10]  was  previously 
calculated  by  Weeks  using  the  CPM.  The  PES  used  is  shown  in  Fig.  5.1.  Like  the  H  +  H2 
problem,  the  reactant  and  product  Mpller  states  had  to  be  calculated.  The  same  procedure 
for  calculating  the  Mpller  states  of  the  H  +  H2  collision  was  used  in  the  Mpller  states 
of  the  OH  +  CO  collision.  These  Mpller  states  are  then  used  to  calculate  the  correlation 
function.  Again,  the  correlation  function  data  was  imported,  and  the  scattering  WDF  was 
calculated  exactly  in  the  same  manner  as  the  H  +  H2  problem.  The  only  difference  was 
that  the  OH  +  CD  correlation  function  had  to  be  down  sampled  while  the  H  +  H2  did 
not.  Because  only  one  arrangement  channel  was  imported,  only  one  scattering  WDF  was 
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calculated.  Calculating  the  scattering  WDF  for  the  OH  +  CO  problem  has  never  been 
computed  previously,  so  the  purpose  was  to  produce  the  OH  +  CO  scattering  WDF.  Due  to 
the  very  sharp  resonant  structures  in  the  transmission  coefficient,  the  OH  +  CO  scattering 
WDF  should  consist  of  unique  features.  The  parameters  and  grid  sizes  associated  with  the 
correlation  functions  are  shown  in  Table  5.1. 


Table  5.1:  H  +  H2  and  OH  +  CO  Parameters 


variable 

H  +  H2 

OH  +  CO 

meaning 

^time 

16000 

1080500 

number  of  time  grid  points 

Hr 

8192 

65536 

number  of  tau  grid  points 

ds 

2 

20 

number  of  points  down  sampled 

^min 

-7999 

-1079999 

minimum  time 

tfnax 

500 

8000 

maximum  time 

The  WDFs  for  the  H  +  H2  and  the  OH  +  CO  collisions  are  different  than  the  square 
potential  problem.  In  the  H+H2  scattering  WDFs,  a  large,  positively  valued,  peak  structure 
still  exists;  however,  this  peak  structure  is  more  localized  than  the  spread  out  chirp  in  the 
square  potentials.  A  series  of  bands  running  nearly  parallel  to  the  time  axis  are  also  present 
in  the  scattering  WDFs.  The  bands  seem  to  diverge  from  the  energy  value  associated  with 
the  large  peak  and  converge  towards  an  asymptotic  upper  energy  limit  or  an  asymptotic 
lower  energy.  These  asymptotic  energies  seem  to  be  centered  on  the  width  of  the  large, 
positive  peak.  The  scattering  WDFs  for  the  upper  two  states  have  a  similar  large  peak 
and  band  structure  that  appears.  The  difference  between  the  ground  state  scattering  WDF 
and  the  upper  two  states  is  that  the  bands  in  the  upper  two  states  converge  asymptotically 
towards  the  energy  value  where  the  peak  value  of  the  scattering  WDF  lies.  Although  the 
ground  state  scattering  WDF  is  approximately  7.5  times  larger  than  the  first  vibrational 
state  scattering  WDF  or  approximately  600  times  larger  than  the  second  vibrational  state. 
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(a)  Scattering  WDF  for  for  v  =  0  ^  v'  =  0 
internal  vibrational  state. 


(b)  Scattering  WDF  for  for  v  =  0  — >  v'  =  1 
internal  vibrational  state. 
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(c)  Scattering  WDF  for  for  v  =  0 
internal  vibrational  state. 
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(d)  Figure  showing  all  three  scaled  scattering 
WDFs  for  V  =  0  ^  v'  =  0,v  =  0  =  1, 

and  V  =  0  ^  v'  =  2  on  the  same  axis 


Figure  5.2:  Scattering  WDFs  for  the  first  three  internal  vibrational  states  of  the  H  +  H2 
collision.  The  reactant  M0ller  state  was  placed  on  the  ground  vibrational  state  and 
correlated  with  product  Mpller  states  on  the  first  three  LSTH  surfaces  labeled  by  the 
internal  vibrational  state,  v.  Each  scattering  WDF  contains  a  large,  positive  valued  single 
peak  and  a  series  of  oscillating  bands  that  seem  to  converge  towards  certain  energy  values. 
A  fourth  plot  contains  all  three  scattering  WDFs  on  the  same  axis.  The  v  =  0  ^  v  =  1 
scattering  WDF  was  scaled  by  a  factor  of  7.95  and  thev  =  0  ^  v  =  2  scattering  WDF  was 
scaled  by  a  factor  of  645.92  such  that  the  peaks  were  equal 


an  interesting  observation  occurs  when  the  three  scattering  WDFs  are  plotted  on  the  same 
axis  and  scaled  such  that  the  large  central  peaks  have  the  same  peak  value.  The  three 
large,  single  peaks  are  above  each  other  on  the  energy  scale.  The  upper  bands  of  the 
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ground  vibrational  state  seem  to  eonverge  with  the  lower  bands  of  the  first  vibrational 
state  to  the  eenter  of  the  first  vibrational  seattering  WDF.  Although  the  upper  bands  of 
the  first  exeited  vibrational  state  eonverged  towards  the  eenter  of  the  peak  of  the  first 
exeited  vibrational  state  and  although  the  lower  bands  of  the  seeond  exeited  vibrational  state 
eonverged  towards  the  eenter  of  the  peak  of  the  seeond  exeited  vibrational  state,  together 
these  bands  form  a  region  very  similar  to  the  region  between  the  ground  state  and  the  first 
vibrational  state.  Further  analysis  is  required  in  order  to  determine  what  these  bands  mean. 
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Figure  5.3:  Seattering  WDF  for  the  OH{v  =  0)  -l-  CO  OCO(v  =  0)  +  H  eollision. 
The  seattering  WDF  eontains  very  long,  narrow  bands  that  eontain  oseillating  peaks  and 
valleys.  The  frequeney  of  oseillations  of  peaks  vary  among  the  different  bands. 


The  OH  +  CO  seattering  WDF  has  a  series  of  very  narrow  band  struetures  that 
run  parallel  to  the  time  grid  for  very  long  periods  of  time  as  eompared  to  the  H  +  H2 
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scattering  WDFs.  There  seems  to  be  a  group  of  bands  that  are  relatively  close  to  each 
other;  however,  between  each  group  of  bands  there  is  a  large  region  of  zero  scattering 
WDF  values.  For  example,  there  is  a  group  of  closely  spaced  bands  at  approximately  0.023 
Hartree  followed  by  another  set  of  closely  spaced  bands  at  approximately  0.0245  Hartree. 
This  is  consistent  with  the  transmission  coefficient  since  the  transmission  coefficient  has  a 
group  of  closely  spaced  sharp  peaks  that  are  spaced  relatively  far  from  the  next  grouping 
of  peaks.  Every  band  in  the  OH  +  CO  data  set  contains  oscillating  positive  and  negative 
peaks;  however,  the  width  of  the  peaks  or  the  magnitude  of  the  peaks  vary.  The  very  long 
bands  arise  from  the  presence  of  long  lived  quasi-bound  states  trapped  by  the  well  of  the 
OH(v  =  0)  +  CO  OCO(v  =  0) -I-//  channel.  Further  analysis  should  be  done  to  determine 
if  the  frequency  of  oscillations  in  the  peaks  of  each  band  (i.e.  the  width  of  each  peak)  have 
any  relation  to  the  parameters  of  the  well  of  OH  +  CO  collision.  Further  analysis  should 
be  done  to  see  if  the  negative  values  of  the  scattering  WDF  correspond  to  correspond  to 
positive  values  in  the  same  location  for  higher  vibrational  states. 

The  H  +  H2  collision  is  an  example  where  a  large  single  peak  is  significantly  much 
larger  than  the  smaller  oscillating  peaks.  As  a  result,  when  the  scattering  WDF  is  projected 
onto  the  frequency  domain,  the  large  single  peak  will  contribute  most  to  |c(w)|^  as  shown 
below  while  the  oscillating  peaks,  which  will  mostly  cancel  each  other  out,  will  have  only  a 
small  contribution  to  |c(w)|^.  Therefore,  the  transmission  coefficient  will  resemble  Fig.  5.4. 
On  the  other  extreme,  the  HO  -v  CO  collision  does  not  have  a  large  peak  but  instead  long 
narrow  bands  of  oscillating  peaks  that  run  parallel  to  the  time  grid.  These  bands  are  the 
major  contributions  |c(w)p.  Hence,  the  transmission  coefficient  for  OH  -l-  CO  consists  of 
many  resonant  sharp  peaks  compared  to  the  H  +  H2  peaks. 
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Figure  5.4:  The  plots  of  the  projection  of  the  scattering  WDF  (blue)  onto  frequency  domain 
as  compared  to  |c(m)^|  respectively  (red)  for  the  H  +  H2{v  =  0)  +  H2{v  =  0)  collision. 


Figure  5.5:  The  plots  of  the  projection  of  the  scattering  WDF  (blue)  onto  frequency  domain 
as  compared  to  |c(m)^|  respectively  (red)  for  the  OH{v  =  0)  +  CO  OCO(v  =  0)  +  H 
collision. 
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VI.  M  +  Ng  Scattering 


With  confidence  that  the  code  is  working  and  with  some  insight  into  the  scattering 
WDF,  we  now  focus  on  producing  scattering  WDFs  for  M  +  Ng  collisions. 

6.1  Theory 

The  Hamiltonian  for  the  M  +  Ng  reaction  is 

1 


H  =  -Pi  +  +  Vf^(rN.,R) 


2fi 


(6.1) 


where  M  is  the  alkali  metal,  Ng  is  the  noble  gas,  yu  is  the  reduced  mass,  7?^  is  the  coordinate 
vector  from  the  origin  to  the  center  of  mass  of  the  two  nuclei,  R  is  the  coordinate  vector 
between  the  center  of  masses  for  each  atom,  is  the  coordinate  vector  from  the  center 
of  mass  of  the  designated  atom  to  the  electron  of  the  designated  atom,  is  the  spin- 
orbit  potential  for  the  alkali  atom  with  quantum  numbers  I  and  s  (the  electronic  angular 
momentum  and  spin  values  respectively),  is  the  momentum  of  the  reduced  mass, 
and  =  HljirNu)  +  ^vg(4^,)  +  VMNgirN^,,  rN;,^,R)  contains  the  electronic 

Hamiltonians  and  the  interaction  potential.  Note  that  the  functional  dependence  of  the  first 
Nm  -  1  electrons  of  alkali  metal  and  the  first  NNg  -  1  electrons  of  noble  gas  were  dropped 
but  still  implied.  The  spin-orbit  potential  for  the  noble  gas  is  zero  because  in  the  ground 
state,  the  total  angular  momentum  for  a  noble  gas  is  zero.  The  total  angular  momentum  is 
zero  because  the  shell  is  full  of  electrons  [22].  The  advantage  of  writing  the  Hamiltonian 
in  the  form  of  Eq.  6. 1  is  that  the  electrostatic  potentials  are  grouped  together  which  makes 
it  easier  to  make  the  Bom-Oppenheimer  expansion. 

Because  of  the  interaction  potential,  the  Hamiltonian  in  Eq.  6.1  is  not  separable 
[4,  5,  23].  In  order  to  resolve  this  problem,  the  Bom-Oppenheimer  approximation  is 
made.  The  Bom-Oppenheimer  approximation  assumes  that  due  to  the  mass  of  the  nuclei 
compared  to  the  mass  of  the  electrons,  the  electrons  will  quickly  settle  into  their  eigenstates 
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every  time  the  nuclei  move.  As  a  consequence,  the  separation  distance,  between  each 
nuclei  is  parametrized.  The  eigenvalues  of  the  electronic  Hamiltonian, 
are  solved  for  each  value  of  R  [20] .  In  other  words,  the  electronic  Hamiltonian  is  completely 
diagonalized,  and  the  diagonal  entries,  the  eigenvalues,  are  plotted  as  a  function  of 
The  plot  of  the  eigenvalue  of  the  Hamiltonian  forms  a  surface  of  potential  energy 
values  called  the  potential  energy  surface  (PES)  [8].  More  specifically,  these  plots  are 
refereed  to  as  the  adiabatic  potential  energy  surfaces  because  the  electronic  Hamiltonian 
is  diagonalized.  A  scattering  event  is  simplified  with  the  Bom-approximation  because  the 
electronic  dynamics  is  eliminated  from  the  problem. 

Unfortunately,  the  Bom-Oppenheimer  approximation  does  not  allow  coupling  be¬ 
tween  the  various  PES.  In  order  to  account  for  this,  a  variation  of  the  Born-Oppenheimer 
approximation,  the  group  Bom-Oppenheimer  approximation,  instead  neglects  non-coupled 
states,  and  the  dimensionality  of  a  coupling  matrix  used  in  the  Born-Oppenheimer  approxi¬ 
mation  reduces  to  the  number  of  coupled  states  with  the  Bom-Oppenheimer  approximation 
being  the  limit  in  which  there  is  a  single  state.  The  premise  for  this  approximation  arises 
from  the  fact  that  the  coupling  terms  can  be  neglected  when  the  adiabatic  electronic  energy 
levels  are  far  apart;  however,  if  the  adiabatic  electronic  energy  levels  approach  each  other  at 
a  particular  nuclear  separation,  then  the  coupling  terms  will  be  large.  Because  this  work  is 
interested  in  the  application  of  scattering  in  a  DPAE  laser  system,  coupling  will  be  assumed 
for  the  group  of  molecular  states  in  the  P-maniforld.  The  molecular  states  of  interest  are: 

3  3  3  3  1  1 

lo’I  In’i  ln’!i  Imi  ln’!i  ^  where  J  is  the  total  angular  momentum  which 

is  the  sum  of  the  total  electronic  angular  momentum  j  and  the  internuclear  angular  momen¬ 
tum  and  is  the  projection  of  J  onto  the  body-fixed  z  axis  [8, 17].  The  body-fixed  states 
are  represented  as  ,4^  >;  however,  the  dependence  on  ^  was  dropped  for  notational  sim¬ 
plicity.  Eor  these  states,  Q  =  +|  when  nij  =  +|  and  Q  =  +|  when  nij  =  ±^.  Because, 
there  are  six  electronic  molecular  states  of  interest,  the  Hamiltonian  will  consist  of  a  set  of 
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6x6  matrices  for  each  value  of  J.  We  were  interested  in  the  body-fixed  projections  rather 
than  the  space-fixed  projections  because  transforming  into  the  body-fixed  reference  frame 
reduces  the  two  body  problem  into  a  one  dimensional  problem  since  the  intemuclear  vector 
^  remains  invariant  [17].  If  we  neglect  the  coupling  between  the  positive  and  negative  z- 
axis  projections  in  the  body-fixed  reference  frame,  the  6x6  matrix  reduces  to  a  3x3  matrix 
where  the  manifold  becomes  two-fold  degenerate. 

The  Hamiltonian  represented  in  the  body-fixed  Born-Oppenheimer  molecular  basis 
[8,  9]  is 
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where  a(R)  is  the  spin-orbit  split  parameter  of  the  alkali  metal  and  where  H  and  S  are 
the  eigenvalues  associated  with  the  adiabatic  electronic  energy  surfaces  labelled  H  and 
S  respectively.  Note  that  the  nij  =  values  are  decoupled  from  the  nij  =  ^ 

values.  With  the  full  Hamiltonian  represented  in  the  Born-Oppenheimer  molecular  basis, 
the  physics  of  the  M  -l-  Ng  can  be  deduced.  The  first  matrix  of  Eq.  6.2  contains  the  radial 
derivative  coupling  terms,  which  are  neglected  [8,  9,  18],  so  the  radial  derivative  coupling 
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matrix  simply  contains  the  kinetic  energy  on  the  diagonal  elements.  The  second  matrix 
in  Eq.  6.2  contains  the  spin-orbit  coupling  terms.  This  matrix  contains  both  the  electronic 
Hamiltonian  and  the  spin-orbit  potential.  The  presence  of  the  spin-orbit  potential  forced 
the  states  with  the  same  electronic  angular  momentum  projected  onto  the  internuclear  axis 
but  originating  from  different  asymptotic  curves  to  couple.  The  spin-orbit  coupling,  which 

3  1 

is  also  known  as  radial  coupling,  couples  the  (the  >  states)  and  the  ^Pi_  (the  > 
states).  The  spin-orbit  coupling  is  the  only  source  of  probability  transfer  at  this  level  of 
theory  between  the  P-manifold  split.  This  coupling  factor  depends  upon  the  strength  of 
the  spin-orbit  coupling  function  and  the  energy  difference  between  the  two  states.  If  the 
energy  of  the  system  is  less  than  the  energy  difference  between  and  ^P^,  then  there  is 
low  probability  that  a  transition  will  occur. 

Once  there  is  a  probability  density  function  on  the  ^Pi  curve  due  to  spin-orbit 
coupling,  Coriolis  coupling  can  occur  on  the  two-fold  degenerate  ^Pi  manifold.  The  matrix 
containing  the  Coriolis  coupling  term  will  be  referred  to  as  the  total  angular  momentum 
matrix  reserving  the  term  Coriolis  coupling  for  the  off  diagonal  matrix  elements.  The 
diagonal  elements  of  this  matrix  has  the  same  functional  form  as  the  effective  potential  in 
the  hydrogen  atom  Hamiltonian,  so  the  diagonal  matrix  elements  will  be  referred  to  as  the 
effective  potential  (or  centrifugal  barrier).  The  total  angular  momentum  matrix  results  from 
the  nuclear  kinetic  energy  begin  transformed  into  the  body-fixed  coordinate  system.  This 
transformation  moves  the  dynamics  of  the  molecule  into  a  rotating  non-inertial  reference 
frame,  and  the  use  of  a  non-inertial  reference  frame  creates  a  fictitious  potential  known 
as  the  centrifugal  potential  along  the  diagonal  of  the  matrix  [6].  This  potential  only 
appears  in  the  body-fixed  coordinate  system  and  is  purely  a  dynamical  effect  dependent 
on  coordinates.  The  effect  of  rotating  into  the  body-fixed  reference  frame  creates  an 
angular  momentum  dependent  term  in  the  Hamiltonian  which  couples  the  total  electronic 
momentum  to  the  nuclear  angular  momentum.  This  coupling  term  measures  the  rate  at 
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which  the  molecular  wave  function  rotates  compared  to  the  nuclear  axis.  If  the  molecular 
wave  function  rotates  at  the  same  rate  as  the  nuclear  axis,  then  the  total  angular  momentum 
is  the  same  as  the  electronic  angular  momentum,  and  no  coupling  occurs.  On  the  other 
hand,  if  the  energies  are  large  enough,  the  nuclear  axis  rotates  at  a  rate  comparable  to  the 
rate  at  which  the  electrons  relax  into  their  adiabatic  states.  The  molecule  is  no  longer  in 
an  adiabatic  state,  and  coupling  between  various  PES  occurs  [8].  Coriolis  coupling  occurs 

3  3 

between  the  I3  >  and  the  i  1  >  states.  It  is  possible  then  that  the  wavefunction  gets  trapped 

2  2 

in  the  manifold  due  to  Coriolis  coupling.  This  prevents  some  of  the  wavefunction  to 
scatter  back  down  to  the  1  manifold.  The  strength  of  the  Coriolis  coupling  depends  upon 
the  velocity  at  which  the  colliding  partners  interact.  As  the  nuclear  velocity  increases,  the 
rate  at  which  the  molecular  axis  rotates  increases.  At  large  enough  energies  the  molecular 
wavefunction  will  not  be  able  to  keep  up  with  the  molecular  axis,  and  a  transition  takes 
place  that  changes  the  electronic  projection  onto  the  molecular  axis. 

The  Hamiltonian  in  Eq.  6.2  is  in  diabatic  form.  In  order  to  have  the  Hamiltonian  in 
adiabatic  form,  the  electronic  terms  need  to  be  diagonalized.  A  transformation  matrix, 
Uso,  is  used  to  diagonalize  the  electronic  terms  [8];  however,  this  transformation  matrix 
must  also  be  applied  to  the  kinetic  operator.  In  this  new  basis,  coupling  terms  are  added  to 
the  kinetic  energy  operator  that  wasn’t  present  in  the  Bom-Oppenheimer  molecular  basis. 
These  new  coupling  terms  are  the  spin-orbit  coupling  physics  in  a  different  representation; 
therefore,  the  spin-orbit  coupling  is  sometimes  referred  to  as  radial  coupling  (not  to  be 
confused  with  radial  derivative  coupling). 

6.2  Results 

Eor  the  M+Ng  systems,  the  correlation  functions  were  previously  calculated  by  Eoper 
[9].  The  Hamiltonian  for  these  systems  is  defined  by  Eq.  6.2.  Because  coupling  between 
the  positive  spin  states  with  the  negative  spin  states  is  small,  the  -  ^  term  in  the  Coriolis 
coupling  component  was  assumed  to  be  zero.  As  a  result,  the  Coriolis  coupling  reduced  to 
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a  3x3  matrix  rather  than  a  6x6  matrix.  The  full  Hamiltonian  was  also  reduced  to  a  3  matrix 


from  eliminating  this  coupling  term.  The  initial  state  was  assumed  to  be  in  the  positive  spin 
state.  The  Mpller  states  were  calculated  using  the  asymptotic  Hamiltonian  and  the  3x3 
full  Hamiltonian.  Because  the  Hamiltonian  is  in  the  diabatic  state,  the  propagation  scheme 
in  the  split  operator  method  had  to  change  slightly.  In  order  to  for  Eq.  3.16  to  apply, 
the  potential  operator  must  act  on  a  wavefunction  in  the  adiabatic  basis  rather  than  the 
diabatic  basis.  Thus  a  transformation  matrix  must  be  applied  to  the  wavefunction.  After  the 
potential  operator  acts  on  the  transformed  wavefunction  and  before  the  Fourier  transform 
is  taken,  the  wavefunction  must  be  transformed  back  into  the  diabatic  representation  using 
the  transpose  conjugate  of  the  transformation  matrix.  The  same  transformation  must  occur 
at  the  second  potential  operator  in  the  split-operator  propagation  scheme. 

Calculation  of  the  Mpller  states  was  done  in  the  same  manner  as  the  H  +  H2  and 
OH  +  CO  problems.  Because  the  full  Hamiltonian  depends  upon  the  total  angular 
momentum,  a  reactant  Mpller  state  was  created  for  each  value  of  J.  As  J  increases,  greater 
amounts  of  energy  was  required  to  reach  the  interaction  potential,  so  the  asymptotic  region 
kept  getting  pushed  back  due  to  the  centrifugal  potential.  The  asymptotic  region  for  all 
regions  was  therefore  defined  according  to  the  start  of  the  asymptotic  region  of  the  highest 
J  value.  For  example,  for  the  K  +  He  system,  the  asymptotic  region  for  J  =  0.5  was 
approximately  20  Bohr  versus  400  Bohr  for  J  =  250.5.  The  highest  J  value  used  was  the 
value  of  J  associated  with  the  particle  being  unable  to  reach  the  interaction  region  for  the  a 
system  corresponding  to  energies  of  0  -  0.0075£’//[9].  For  each  system,  the  highest  J-value 
was  different,  so  each  system  had  a  different  number  of  correlation  functions.  Because 

3  3  1 

there  are  three  states  of  interest  for  each  value  of  J:  \{,\  >,  \\,\  >,  and  \\,\  >,  three  product 

2  2  2  2  2  2 

Mpller  states  for  each  value  of  J  were  created  on  each  surface.  Not  only  was  there  a  static 
product  Mpller  state  for  each  surface,  the  reactant  Mpller  states  were  also  propagated  along 
each  curve  and  correlated  with  each  product  Mpller  state.  In  other  words,  for  each  alkali 
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plus  noble  gas  system,  there  were  eorrelation  funetions  eomputed  for  eaeh  of  the  reaetant 
Mpller  states  evolving  on  the  three  different  PES  overlapping  with  eaeh  of  the  three  produet 
Mpller  states  as  a  funetion  of  total  angular  momentum. 

Fig.  6. 1  sumarizes  the  eomputational  proeedure  along  with  the  eoupling  terms  present. 
Although  all  of  the  states  belonging  to  the  maniforld  are  at  the  same  energy  level  in  the 
asymptotie  regime,  they  were  artifieally  seperated  in  Fig.  6.1  in  order  to  show  whieh  states 
leads  to  seperate  adiabatie  surfaees  duing  the  eollision.  The  adiabatie  surfaees  are  shown 
in  Fig.  6.2. 
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Figure  6.1:  Summary  of  the  M  +  Ng  collision.  The  states  13,3  >  and  >  belong  to  the 

2  2  2  2  ^ 

1 

manifold  in  the  asymptotic  regime  while  \{,\>  belongs  to  the  ^Pi  manifold.  Although  all 

22  2 

of  the  states  belonging  to  the  ^Pi  maniforld  are  at  the  same  energy  level  in  the  asymptotic 
regime,  they  were  artifieally  seperated.  Three  product  Mpller  states  are  placed  on  each 

of  the  adiabatic  PES  which  are  labeled  as  Hi  ,  Tls,  and  Zi.  A  reactant  Mpller  can  be 

2  2  2 

placed  on  any  adiabatic  surface,  but  as  an  example,  the  reactant  Mpller  state  starts  on  the 
Tti  surface.  Radial  coupling  couples  the  Hi  and  Zi  surfaces  when  the  reactant  Mpller 
state  is  propagated  into  the  interaction  region,  and  Coriolis  coupling  couples  the  Tti  and  Zi 
surfaces. 


Scattering  WDFs  were  produced  for  the  following  systems:  Ke+He,  Ke+Ne,  K+Ar, 
Rb+He,  Rb+Ne,  Rb+Ar,  Cs+Ne,  Cs+Ar.  The  same  code  that  imported  the  correlation 
functions  and  produced  the  scattering  WDFs  for  previous  potentials  was  also  used  for  the 
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Figure  6.2:  The  adiabatic  PES  (11^  (blue),  IIs  (red),  and  (green)  for  the  collision. 


M  +  Ng  collisions.  After  the  scattering  WDFs  were  calculated,  an  analysis  of  the  scattering 
WDFs  was  done.  The  parameters  and  grid  sizes  associated  with  the  correlation  function  of 
each  alkali  plus  noble  gas  studied  are  shown  in  Table  6.1. 


Table  6.1:  Alkali  and  Noble  Gas  Parameters 


meaning 

K  + He 

K  +  Ne 

K  +  Ar 

Rb  + He 

Rb  +  Ne 

Rb  +  Ar 

Cs  +  Ne 

Cs+Ar 

time  grid 

150000 

125000 

125000 

100000 

75000 

125000 

15000 

125000 

tau  grid 

65536 

65536 

65536 

65536 

65536 

65536 

65536 

65536 

sampled 

2 

2 

2 

2 

1 

2 

1 

2 

min  time 

0 

0 

0 

0 

0 

0 

0 

0 

max  time 

2999980 

4999960  4999960  1999980 

5999920 

4999960 

5999920 

4999960 

Figs.  6. 3-6. 5  show  the  scattering  WDFs  for  the  K  -l-  Ne  system  for  the  case  when 

1 

the  reactant  Mpller  state  was  on  the  \]  >  surface  and  correlated  with  product  states  on 

2 

each  surface  at  a  total  angular  momentum  value  of  7  =  100.5.  As  we  will  see  shortly. 
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the  ’’reflection  scattering  WDF”,  which  is  the  scattering  WDF  associated  with  the  product 

Mpller  state  on  the  same  surface  as  the  reactant  Mpller  state,  had  the  greatest  amplitude 

and  maintained  the  same  structure  as  the  total  angular  momentum  increased.  The  structure 

of  this  scattering  WDF  consisted  of  closely  spaced,  negatively  chirped,  single  peaked  (or 

valley)  bands.  The  bands  alternate  between  positive  and  negative  values.  The  scattering 

1 

WDF  for  the  correlation  function  between  a  reactant  Mpller  state  on  the  |  ^  >  surface  and 

2 

3  3 

product  Mpller  state  on  the  |(  >  or  the  I3  >  surface  also  have  a  series  of  closely  space 

2  2 

bands  with  a  negative  chirp  to  their  structure;  however,  unlike  the  scattering  WDF  for 

1  I 

the  Ci’i(t),  some  of  the  bands  consist  of  oscillating  peaks  and  valleys  inside  them  rather 

2’2 

than  being  a  single  peak.  If  the  square  well  potential  provides  insight  into  the  meaning 
of  these  oscillating  peaks  and  valleys  that  appear  in  the  middle  of  the  scattering  WDF  for 
the  M  +  Ng  collisions,  the  appearance  of  the  oscillating  peaks  and  valleys  in  the  band  may 
signify  multiple  transmissions  of  probability  at  that  surface  depending  upon  the  number  of 
bands  present  with  oscillating  peaks  of  valleys.  The  reflection  scattering  WDF  has  a  similar 
structure  to  the  square  well  scattering  WDF  at  low  square  well  depths  which  had  only 
a  single  transmission  of  probability  through  the  square  well  with  significant  probability. 
The  bands  with  oscillating  peaks  and  valleys  did  not  appear  until  the  square  well  depth 
was  deep  enough  such  that  any  secondary  transmissions  (or  reflections)  had  significant 
probability.  This  may  apply  that  there  were  multiple  transfers  of  relatively  significant 

3  3 

probability  amplitude  to  the  1^  >  or  the  I3  >  surfaces.  As  the  total  angular  momentum 

2  2 

changed,  the  structure  of  the  scattering  WDF  may  change;  but  the  scattering  WDF  still  had 
these  basic  features. 


6.2.1  Detailed  Balance. 

When  the  total  angular  momentum  is  equal  to  1.5  atomic  units,  J  =  1.5,  the  Coriolis 
coupling  and  effective  potential  are  very  small.  The  only  coupling  that  occurs  is  the  radial 
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Figure  6.3:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant  M0ller 
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state  in  1 1  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.4:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant  M0ller 
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state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.5:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant  M0ller 
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state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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coupling,  and  if  the  reactant  Mpller  state  starts  on  either  the  >  or  the  \  ]  >  PES,  then  the 

2  2 

dynamics  of  probability  transfer  can  be  estimated  by  a  two  level  system  instead  of  three. 

According  to  the  principle  of  detailed  balance,  the  rate  of  probability  transfer  between  the 

1 

two  surfaces  is  the  same  when  the  reactant  Mpller  state  starts  on  either  the  \]  >  or  the 

2 

3 

1 1  >  PES  up  to  a  proportionality  constant  governed  by  the  difference  in  energy  levels.  The 

2 

principle  of  detailed  balance  is  shown  below: 


Pij(E) 

Pji(E) 


_  AE 

=  e  et 


(6.3) 


1  3 

where  AE  is  the  difference  in  energy  levels  between  the  \]  >  and  the  li  >  surfaces,  K  is 

2  2 

the  Boltzmann  constant,  T  is  the  temperature,  and  is  the  integrated  projection  of  the 


scattering  WDE  at  energy  E  with  the  reactant  Mpller  state  starting  on  surface  i  and  ending 


on  surface  j.  Eigs.  6.6-6.11  demonstrate  this  idea.  Notice  that  no  detail  of  the  WDE  can 

3  13 

be  seen  on  the  I3  >  in  this  scale,  but  the  WDE  can  be  seen  on  both  the  II  >  and  the  II  > 

2  2  2 

1 

surface.  In  Eig.  6.6,  the  reactant  Mpller  state  started  on  the  1 1  >  surface  of  .ST  -1-  Ne,  so  most 


of  the  probability  reflected  back  on  the  |  j  >  surface  while  about  a  tenth  of  the  probability 

2 

3 

was  coupled  via  radial  coupling  to  the  |  ]  >  surface  as  seen  by  the  amplitudes  of  the  WDE 

2 

for  each  surface.  Now  Eig.  6.7  shows  the  WDEs  for  the  K  -i-  Ne  system  when  the  reactant 

3 

Mpller  state  starts  on  the  \]  >  surface.  As  the  figure  shows,  most  of  the  probability  is 

2 

3 

reflected  back  on  the  |  j  >  surface  while  about  a  tenth  of  the  probability  is  transferred  to  the 

2 

1 

1  (  >  surface  via  radial  coupling. 

2 

1 

When  the  WDE  associated  with  |i  >  product  Mpller  state  in  Eig.  6.6  is  compared  with 

2 

3 

the  Ij  >  surface  in  Eig.  6.7,  there  is  a  lot  of  similarity  between  the  scattering  WDEs.  The 

2 

similarity  between  Eig  6.8  and  Eig  6.9  is  consistent  with  the  principle  of  detailed  balance. 


Since  Pij(E)  differes  from  Pji(E),  this  suggests  that  the  scattering  WDEs  differ  by  a  constant 


amplitude  factor,  e  kt  .  Eurther  analysis  must  be  done  in  order  to  determine  if  the  difference 


in  amplitude  is  e  .  When  the  WDE  associated  with  II  >  product  Mpller  state  in  Eig.  6.8 


is  compared  with  the  \]  >  surface  in  Eig.  6.9,  the  similarities  are  also  consistent  with  the 

2 
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3  3 

principle  of  detailed  balance.  The  Cs+Ne  eollision  between  the  |  j  >  and  1 1  >  surfaee  is  also 

2  2 

eonsistent  with  the  principle  of  detailed  balance  sinee  the  large  spin-obit  split  of  Ci’  -l-  Ne 


eauses  the  radial  eoupling  to  be  insignifieant  ereating  a  two-level  system.  An  example  of 


this  is  shown  for  Cs  +  Ne  at  J  =  150.5  in  Fig.  6.12-Fig.  6.17.  Figs.  6.6-6.17  suggest 


that  the  principle  of  detailed  balanee  is  not  only  satisfied  for  the  integrated  projeetion  of 


the  seattering  WDF  at  energy  E  for  a  two  level  system  but  throughout  the  entire  collision. 


Again  further  analysis  must  be  done  to  verify  this  claim. 


6.2.2  High  Total  Angular  Momentum  Values. 


At  high  values  of  total  angular  Momentum,  the  effective  potential  becomes  so  large 


that  the  radial  (spin-orbit)  eoupling  is  essentially  turned  off  as  seen  in  Eq.  6.2.  The  effeetive 


potential,  which  are  the  diagonal  elements  of  the  angular  momentum  dependent  matrix  in 


Eq.  6.2,  shifts  the  adiabatie  PES  to  higher  values  of  internuelear  separation  distances,  R. 


Beeause  the  spin-orbit  eoupling  function  is  independent  of  total  angular  momentum,  the 


spin-orbit  coupling  function  remains  the  same  which  radially  couples  the  |  i  >  and  the 

2 

3 

Ij  >  surfaces  localized  at  small  values  of  internuelear  separation  distances  (between  5-10 

2 

Bohr  approximately).  When  the  effeetive  potential  shifts  the  PES  to  higher  values  of  R, 


only  high  energy  eomponents  of  the  wave-packet  will  penetrate  deep  enough  into  the  PES 


and  experience  the  spin-orbit  eoupling  effeets  while  the  low  energy  eomponents  will  never 


experienee  this  coupling  effect. 


If  the  reactant  Mpller  state  starts  on  the  |  ]  >  surface,  then  the  seattering  WDE  should 

2 

3  3 

decrease  in  amplitude  for  the  |  I  >  and  the  I3  >  surfaees.  The  radial  Coupling  gets  turned 

2  2 

off  for  most  values  of  energy  being  studied,  so  most  of  the  probability  is  reflected  back  on 

1 

the  I  ]  >  surfaee  before  the  effeets  from  spin-orbit  eoupling  ean  be  felt.  Beeause  very  little 

2 

3 

probability  is  transferred  to  the  i  (  >  surfaee,  there  is  very  little  probability  to  be  transferred 


to  the  1 1  >  surfaee  via  Coriolis  eoupling.  Eigs.  6.18-6.24  show  that  as  J  increases,  the 

2 


82 


amplitude  of  the  scattering  WDF  decreases  for  the  |  1  >  and  the  I3  >  surfaces  when  the 

2  2 

1 

reactant  Mpller  state  starts  on  the  |  J  >  surface  of  ^  +  Ne. 

2 

3 

The  scale  used  in  Fig.  6.23  shows  that  the  scattering  WDF  is  zero  for  the  |  1  >  and 

2 

3 

the  1|  >  surfaces;  but,  when  the  scale  is  changed  to  0.3  which  is  approximately  -j^  of  the 

2 

3 

original  scattering  WDF  amplitude  on  the  |i  >  surface  when  the  effective  potential  and 

2 

Coriolis  coupling  were  turned  off  (J  =  1.5),  the  scattering  WDF  can  be  seen  as  shown 
in  Fig.  6.24.  Also,  one  can  observe  the  scattering  WDF  shift  to  higher  energy  levels  as  J 

3  3 

increases  for  the  1 1  >  and  the  I3  >  surfaces  in  Fig.  6.25-Fig.  6.30. 

2  2 

3  3 

If  the  reactant  Mpller  state  starts  on  either  the  1 1  >  or  the  1 3  >  surface,  the  scattering 

2  2 

1 

WDF  should  decrease  in  amplitude  for  the  |  ^  >  surface  for  the  same  reasons  mentioned 

2 

above  as  shown  in  Figs.  6.31-6.36.  The  radial  coupling  term  is  turned  off,  and  the 

3  3 

probability  is  only  transferred  between  the  |  I  >  and  the  1 3  >  surfaces  (the  Coriolis  coupled 

2  2 

surfaces)  as  shown  in  Figs.  6.37-6.48.  At  low  values  of  J,  the  Coriolis  coupling  is  weak,  so 

3  3 

there  isn’t  much  probability  that  flows  between  the  1 1  >  and  the  I3  >  surfaces  regardless  of 

2  2 

where  the  reactant  Mpller  state  starts.  As  J  increases,  the  strength  of  the  Coriolis  coupling 

3 

increases  which  increases  the  amplitude  of  the  scattering  WDF  on  the  \]  >  surface  if  the 

2 

3  3 

reactant  Mpller  state  started  on  the  1 1  >  surface  or  the  1 3  >  surface  if  the  reactant  Mpller 


state  started  on  the  |  ]  >  surface.  For  high  values  of  J,  radial  coupling  is  turned  off  at  large 

2 

3  3 

J,  so  most  of  the  probability  (all  but  a  fraction  of  a  percent)  exists  on  the  |  ?  >  and  the  1  !  > 

2  2 

surfaces.  The  amplitude  of  the  scattering  WDF  for  these  surfaces  are  on  the  same  order  of 
magnitude  for  high  values  of  J.  Similarly,  the  structure  of  the  scattering  WDFs  are  almost 

3  3 

the  same  for  high  values  of  J  for  the  |  I  >  and  the  I3  >  surfaces.  This  suggests  that  the 

2  2 

3  3 

probability  transfer  between  the  1 1  >  and  the  1 1  >  surfaces  reaches  some  sort  of  equilibrium 

2  2 

for  large  values  of  J. 
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Figure  6.6:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant  M0ller 

1  1 

state  in  1 1  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.7:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant  M0ller 

3  3 

state  in  1 1  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.8:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant  M0ller 

1  3 

state  in  1 1  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.9:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant  M0ller 

3  1 

state  in  1 1  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.10:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |  1  >’  product  Mpller  state  in  I3  >,  and  J  =  1.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.8.  The  purpose  of  this  is  to  show  that  collision  reduces  to  a  two  level 
system  at  low  values  of  J. 
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Figure  6.11:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  1  >’  product  Mpller  state  in  I3  >,  and  J  =  1.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.8.  The  purpose  of  this  is  to  show  that  collision  reduces  to  a  two  level 
system  at  low  values  of  J. 
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Figure  6.12:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  150.5 
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Figure  6.13:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  3  >’  product  Mpller  state  in  I3  >,  and  J  =  150.5 
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Figure  6.14:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  i  >’  product  Mpller  state  in  I3  >,  and  J  =  150.5 
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Figure  6.15:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  3  >’  product  Mpller  state  in  | ,  >,  and  J  =  150.5 


2 


2 


93 


(riT) 


o  o  o  o 

oooo  oooo 

OOOO  f-C\JCO'^ 

CO  C\J  ■»-  O  I  I  I  I 


Figure  6.16:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  1 

M0ller  state  in  |  1  >’  product  Mpller  state  in  1]  >,  and  J  =  150.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.1 12.  The  purpose  of  this  is  to  show  that  collision  reduces  to  a  two  level 
system  at  low  values  of  J. 
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Figure  6.17:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  1 

M0ller  state  in  ||  >,  product  Mpller  state  in  1]  >,  and  J  =  150.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.1 12.  The  purpose  of  this  is  to  show  that  collision  reduces  to  a  two  level 
system  at  low  values  of  J. 
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Figure  6.18:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.19:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |  1  >’  product  M0ller  state  in  |  i  >’  and  J  =  50.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.18.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.20:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  >,  product  Mpller  state  in  |i  >,  and  J  =  100.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.18.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.21:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |  1  >’  product  Mpller  state  in  \]  >,  and  J  =  150.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.18.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 


99 


Time  (A.U) 


Ln 


O 

X 


Figure  6.22:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  >,  product  Mpller  state  in  |i  >,  and  J  =  200.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.18.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.23:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  >,  product  Mpller  state  in  |i  >,  and  J  =  250.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.18.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.24:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  |  i  >’  and  J  =  250.5.  The  scale  for  this 

2  2 

figure  is  different  than  the  scale  used  in  Fig.  6.18.  The  scale  is  based  on  the  maximum  and 
minimum  value  of  the  scattering  WDF  to  show  that  the  amplitude  exists. 
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Figure  6.25:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.26:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |  1  >’  product  Mpller  state  in  I?  >,  and  J  =  50.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.25.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.27:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  >,  product  Mpller  state  in  I3  >,  and  J  =  100.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.25.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.28:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |  1  >’  product  Mpller  state  in  I3  >,  and  J  =  150.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.25.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.29:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  >,  product  Mpller  state  in  I3  >,  and  J  =  200.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.25.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.30:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  >,  product  Mpller  state  in  I3  >,  and  J  =  250.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.25.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.31:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  1 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  0.5 
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Figure  6.32:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  1 

M0ller  state  in  |  1  >’  product  M0ller  state  in  |  i  >’  and  J  =  50.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.31.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.33:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  1 

M0ller  state  in  >,  product  Mpller  state  in  |i  >,  and  J  =  100.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.31.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.34:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  1 

M0ller  state  in  |  1  >’  product  Mpller  state  in  \]  >,  and  J  =  150.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.31.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.35:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  1 

M0ller  state  in  >,  product  Mpller  state  in  |i  >,  and  J  =  200.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.31.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.36:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  1 

M0ller  state  in  >,  product  Mpller  state  in  |i  >,  and  J  =  250.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.31.  The  purpose  of  this  is  to  observe  the  diminishing  amplitude  of  the 
scattering  WDF  as  the  J  increases. 
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Figure  6.37:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  i  >’  product  Mpller  state  in  ||  >,  and  J  =  0.5 
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Figure  6.38:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  50.5 
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Figure  6.39:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  | ,  >,  and  7  =  100.5 
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Figure  6.40:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  150.5 
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Figure  6.41:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  200.5 
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Figure  6.42:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  250.5 
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Figure  6.43:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  1  >’  product  Mpller  state  in  I3  >,  and  J  =  0.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.37.  The  purpose  of  this  is  to  observe  the  increasing  amplitude  of  the 
scattering  WDF  as  J  increases  and  to  show  how  the  amplitude  of  the  scattering  WDFs  for 

3  3 

the  I!  >  and  the  I3  >  surfaces  approach  each  other  for  high  J. 
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Figure  6.44:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  1  >’  product  Mpller  state  in  I?  >,  and  J  =  50.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.37.  The  purpose  of  this  is  to  observe  the  increasing  amplitude  of  the 
scattering  WDF  as  J  increases  and  to  show  how  the  amplitude  of  the  scattering  WDFs  for 

3  3 

the  I!  >  and  the  I3  >  surfaces  approach  each  other  for  high  J. 

2  2 


122 


Time  (A.U) 


Figure  6.45:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  >,  product  Mpller  state  in  I3  >,  and  J  =  100.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.37.  The  purpose  of  this  is  to  observe  the  increasing  amplitude  of  the 
scattering  WDF  as  J  increases  and  to  show  how  the  amplitude  of  the  scattering  WDFs  for 

3  3 

the  I!  >  and  the  I3  >  surfaces  approach  each  other  for  high  J. 

2  2 
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Figure  6.46:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  1  >’  product  Mpller  state  in  I3  >,  and  J  =  150.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.37.  The  purpose  of  this  is  to  observe  the  increasing  amplitude  of  the 
scattering  WDF  as  J  increases  and  to  show  how  the  amplitude  of  the  scattering  WDFs  for 

3  3 

the  I!  >  and  the  I3  >  surfaces  approach  each  other  for  high  J. 
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Figure  6.47:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  >,  product  Mpller  state  in  I3  >,  and  J  =  200.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.37.  The  purpose  of  this  is  to  observe  the  increasing  amplitude  of  the 
scattering  WDF  as  J  increases  and  to  show  how  the  amplitude  of  the  scattering  WDFs  for 

3  3 

the  I!  >  and  the  I3  >  surfaces  approach  each  other  for  high  J. 
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Figure  6.48:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

3  3 

M0ller  state  in  >,  product  Mpller  state  in  I3  >,  and  J  =  250.5.  Note  that  the  scale  is 

2  2 

the  same  as  Fig.  6.37.  The  purpose  of  this  is  to  observe  the  increasing  amplitude  of  the 
scattering  WDF  as  J  increases  and  to  show  how  the  amplitude  of  the  scattering  WDFs  for 

3  3 

the  I!  >  and  the  I3  >  surfaces  approach  each  other  for  high  J. 

2  2 
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6.2.3  Angular  Momentum  Study. 

Scattering  WDFs  were  produced  for  the  Potassium  plus  noble  gas  collisions  in  steps 


of  A7  =  5.0  starting  from  7  =  5.5  to  7  =  250.5  (7  =  1.5  is  also  included)  with  the 

\ 

reactant  state  starting  in  the  \]  >  surface.  The  purpose  of  this  study  was  to  see  how  the 

2 

total  angular  momentum  impacted  probability  transfer  between  the  various  surfaces.  As 
mentioned  in  the  square  well  problem,  the  sum  of  all  the  values  of  the  scattering  WDF 
along  an  energy  slice  multiplied  by  a  normalization  factor  is  the  probability  of  transmitting 
on  that  surface.  The  normalization  factor  is  shown  in  Eq.  2.10  and  is  independent  of  the 
value  of  7.  Although  the  maximum  and  minimum  values  of  the  scattering  WDF  along  the 
energy  slice  are  only  two  contributions  to  the  sum,  their  difference  provides  the  largest 
contribution  to  the  sum.  As  discussed  in  the  square  well  problems,  for  a  given  peak 
there  is  an  associated  valley  on  the  same  order  of  magnitude  with  the  same  peak  width, 
so  they  cancel  each  other  out  in  the  sum.  The  largest  peak  and  the  largest  valley  have 
the  greatest  difference  in  their  magnitude  and  do  not  necessarily  cancel  each  other  out,  so 
the  difference  between  the  largest  peak  and  valley  can  be  used  to  approximate  the  sum 
of  the  scattering  WDF  values  along  a  particular  energy  slice.  Furthermore,  the  maximum 
and  minimum  values  for  a  given  energy  slice  are  approximately  the  same  for  a  range  of 
energy  values  (multiple  energy  slices),  so  we  can  use  the  maximum  and  minimum  values 
of  the  scattering  WDF  to  compare  the  relative  flow  of  probabilities  at  various  total  angular 
momentum  among  the  various  potential  energy  surfaces.  Because  we  are  concerned  with 
the  change  in  probability  transfer  as  7  changes,  we  simply  looked  at  the  maximum  value  of 
the  scattering  WDF  instead  of  the  difference  between  the  maximum  and  minimum  value  of 
the  scattering  WDF  as  shown  in  Fig.  6.49. 

1 

When  the  reactant  Mpller  state  starts  on  the  |  (  >  surface,  there  is  very  little  probability 

2 

3 

on  the  1 3  >  for  low  values  of  7,  and  the  maximum  value  of  the  scattering  WDF  associated 

2 

3 

with  the  product  Mpller  state  on  the  1 1  >  surface  is  very  near  zero  for  low  values  of  7.  As 

2 
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the  Coriolis  coupling  becomes  stronger  as  J  increases,  the  maximum  value  of  the  scattering 

3 

WDF  for  the  I3  >  surface  increases.  The  maximum  value  of  the  scattering  WDF  peaks  and 

2 

decreases  for  even  higher  J  until  the  maximum  value  for  the  scattering  WDF  becomes  zero 
for  very  high  J.  The  maximum  value  of  the  scattering  WDF  peaks  and  decreases  to  zero 
due  to  the  effects  of  the  effective  potential  which  reduces  the  amount  of  probability  that 
can  be  radially  coupled.  The  maximum  value  of  the  scattering  WDF  associated  with  the 

3 

product  Mpller  state  on  the  |  J  >  surface  is  more  complicated  due  to  the  interplay  between 

2 

Coriolis  and  radial  coupling  at  that  surface.  When  the  Coriolis  coupling  is  turned  off  at 

3 

low  values  of  J,  radial  coupling  is  the  only  source  of  probability  to  the  |  ]  >  surface,  so 

2 

3 

unlike  the  1 1  >  surface  scattering  WDF,  the  maximum  value  of  the  scattering  WDF  is  non- 

2 

zero  at  low  values  of  J.  As  J  increases,  the  strength  of  the  Coriolis  coupling  increases 

3 

and  some  of  the  probability  is  transferred  to  the  I3  >  surface.  At  7  =  60.5,  the  effective 

2 

3 

potential  overpowers  the  Coriolis  coupling  (since  the  probability  transfer  to  the  |  ]  >  surface 

2 

decreases  with  increasing  effective  potential),  and  the  maximum  value  of  the  scattering 

3 

WDF  associated  with  the  product  Mpller  state  on  the  >  surface  drops  to  zero  at  the  same 


Figure  6.49:  The  maximum  value  of  the  scattering  Wigner  Distribution  Function  for  K-l-Ne 

1  3 

system  with  reactant  Mpller  state  in  |  I  >’  product  Mpller  state  in  I?  >  (Blue),  and  product 

2  2 
3 

Mpller  state  I3  >  (Red).  The  maximum  value  was  sampled  for  every  AJ  =  5  ranging  from 
7  =  5.5  to  7  =  250.5  (including  7  =  1.5. 
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Figure  6.50:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ]]  >,  and  J  =  10.5 
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Figure  6.51:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  30.5 
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Figure  6.52:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  35.5 
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Figure  6.53:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  60.5 
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Figure  6.54:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  85.5 
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Figure  6.55:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |  j  >,  product  Mpller  state  in  ||  >,  and  7  =  1 15.5 
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Figure  6.56:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  145.5 
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Figure  6.57:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  170.5 
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Figure  6.58:  Scattering  Wigner  Distribution  Function  for  K+Ne  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ]]  >,  and  J  =  225.5 
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An  interesting  feature  that  occurs  on  the  |  ]  >  surface  in  Fig.  6.49  is  the  appearance  of 

2 

a  sharp  dip  in  the  maximum  value  of  the  scattering  WDF  at  7  =  30.5.  This  sharp  dip  is 
followed  by  one  smaller  dip  at  7  =  50.5.  These  dips  do  not  correspond  to  a  peak  on  the 

3 

1 3  >  surface  in  Fig.  6.49;  but,  this  dip  does  follow  a  similar  occurrence  in  the  change  in 

2 

the  structure  of  the  scattering  WDF.  As  shown  in  Figs.  6.50-6.58,  the  number  of  oscillating 
peaks  and  valleys  that  appear  in  the  middle  of  the  structure  of  the  scattering  WDF  increases 
as  7  increases.  When  the  effective  potential  begins  to  turn  off  the  radial  coupling,  the 
number  of  oscillating  peaks  and  valleys  that  appear  in  the  middle  of  the  scattering  WDF 
begins  to  decrease.  By  7  =  225.5,  there  are  no  more  oscillating  peaks  and  valleys  just 
like  at  low  values  of  7.  After  the  dip  found  in  Fig.  6.49,  the  number  of  oscillating  peaks 
and  valleys  in  the  scattering  WDF  decrease.  At  7  =  35.5  which  is  located  in  the  dip,  the 
scattering  WDF  seems  to  have  split  in  half  as  shown  in  Fig.  6.52.  Further  investigation  of 
the  oscillating  peaks  and  valleys  in  the  scattering  WDF  reveals  that  not  only  do  the  number 
of  oscillating  peaks  and  valleys  increase  (decrease),  the  number  of  bands  of  oscillating 
peaks  and  valleys  increase  (decrease)  perhaps  suggesting  multiple  transfers  of  probability 
between  surfaces. 

6.2.4  Noble  Gases. 

If  the  noble  gas  collision  partner  is  changed  while  the  alkali  metal  remains  constant, 
the  angular  momentum  matrix  in  Eq.  6.2  will  be  affected.  The  radial  coupling  will  be 
unaffected  because  the  noble  gas  does  not  contribute  to  spin-orbit  potential.  As  the  mass 
of  the  noble  gas  increases,  the  reduced  mass  increases,  and  the  effective  potential  and  the 
Coriolis  coupling  terms  in  the  angular  momentum  matrix  decrease  in  magnitude  for  a  given 
value  of  7.  As  a  result,  higher  7  values  are  needed  in  order  to  achieve  the  same  Coriolis 
coupling  and  effective  potential  strengths  seen  in  the  lighter  noble  gas  partners  such  as 
Helium.  Figs.  6.59-6.70  demonstrate  this  concept.  The  amplitude  of  the  scattering  WDF 

3 

for  K  +  He  drops  a  couple  orders  of  magnitude  by  the  time  7  =  105.5  for  the  |  I  >  and  the 
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3  1 

1 3  >  surface  if  the  reaetant  M0ller  state  starts  on  the  |  j  >  surfaee;  however,  the  amplitude 

2  2 

of  the  scattering  WDF  for  K  +  Arh  on  the  same  order  of  magnitude  for  the  product  Mpller 

3 

state  on  the  \\  >  surfaee  for  J  =  105.5  as  eompared  to  the  seattering  WDF  amplitude  for 

2 

3 

7  =  1.5  The  amplitude  of  the  seattering  WDF  for  the  product  Mpller  state  on  the  i  3  > 

2 

surfaee  peaked  at  around  J  =  105.5  as  compared  to  the  amplitude  of  the  seattering  WDF 
at  7  =  1.5.  In  fact,  the  amplitude  of  the  seattering  WDF  for  these  surfaees  for  the  K  +  Ar 
eollision  drops  a  couple  orders  of  magnitude  when  7  =  275.5  which  is  approximately  3 
times  the  total  angular  momentum  value  than  the  K  +  He  system. 

Figs.  6.71-6.73  are  the  plots  of  the  maximum  value  of  the  scattering  WDF  as  a  function 

of  7  for  the  K  +  He,  K  +  Ne,  and  K  +  Ar  systems  respeetively  when  the  reactant  Mpller 

1 

state  starts  on  the  |  j  >  surfaee.  As  the  noble  gas  inereases  in  mass,  the  differenee  in  the 

2 

3  3 

order  of  magnitudes  of  the  i  \  >  and  1 1  >  surfaees  inereases.  For  example,  for  K  +  He  the 

2  2 

3 

maximum  value  of  the  seattering  WDF  for  the  1 3  >  is  comparable  to  the  maximum  value  of 

2 

3 

the  1 3  >  WDF;  however,  for  K  +  Ne  and  K  +  Ar  this  is  never  true.  The  maximum  value  of 

2 

3 

the  seattering  WDF  for  the  1 3  >  surface  never  reaches  the  same  order  of  magnitude  as  the 

2 

3 

1 1  >  surfaee  for  K  +  Ne.  The  difference  between  the  relative  order  of  magnitudes  between 

2 

these  two  surfaces  is  slightly  greater  for  the  K  +  Ar  system.  The  figures  also  reveal  that 
the  dip  found  in  Fig.  6.49  for  the  study  of  the  peak  amplitude  of  the  seattering  WDF  as  a 
function  of  7  increases  in  depth  as  the  noble  gas  becomes  heavier. 

3  3 

When  the  reactant  Mpller  state  starts  on  the  |  J  >  or  the  1 1  >  surfaces  for  these  noble 

2  2 

gas  eollisions,  radial  eoupling  is  turned  off  at  high  values  of  7  due  to  the  effeetive  potential 

as  shown  in  Figs.  6.74-6.89.  The  amplitude  of  the  seattering  WDF  for  the  product  state  on 

1  1 

the  1 1  >  surface  deeays  to  zero,  and  there  is  no  probability  transfer  to  the  |  ]  >  state.  As  the 

2  2 

mass  of  the  noble  gas  increases  the  7  value  at  which  radial  coupling  is  turned  off  changes. 
This  implies  that  the  amplitude  of  the  seattering  WDF  for  this  surface  remains  signifieant 
for  larger  values  of  7. 
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Figure  6.59:  Scattering  Wigner  Distribution  Function  for  K+He  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.60:  Scattering  Wigner  Distribution  Function  for  K+He  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  50.5 
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Figure  6.61:  Scattering  Wigner  Distribution  Function  for  K+He  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.62:  Scattering  Wigner  Distribution  Function  for  K+He  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.63:  Scattering  Wigner  Distribution  Function  for  K+He  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  50.5 
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Figure  6.64:  Scattering  Wigner  Distribution  Function  for  K+He  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.65:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.66:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.67:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  200.5 
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Figure  6.68:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.69:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  100.5 
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Figure  6.70:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

1  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  200.5 
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Figure  6.11:  The  maximum  value  of  the  seattering  Wigner  Distribution  Funetion  for  K+He 

1  3 

system  with  reaetant  M0ller  state  in  li  >,  produet  Mpller  state  in  li  >  (Blue),  and  produet 

2  2 
3 

Mpller  state  I3  >  (Red).  The  maximum  value  was  sampled  for  every  AJ  =  5  ranging  from 
7  =  5.5  to  7  =  250.5  (ineluding  7  =  1.5). 


3 

If  the  reaetant  Mpller  state  starts  on  the  1 1  >  surfaee,  there  has  to  be  Coriolis  eoupling 

2 

3  1 

to  the  I  ]  >  surfaee  before  there  ean  be  radial  eoupling  to  the  1  ]  >  surfaee.  By  the  7  value  that 

2  2 

Coriolis  eoupling  is  high  enough  to  allow  for  signifieant  probability  to  be  radially  eoupled 


to  the  I  (  >,  the  effeetive  potential  begins  to  reduee  the  amount  of  probability  going  to  the 

2 

1 

|(  >.  As  shown  in  Figs.  6.87-6.89,  the  amplitude  of  the  seattering  WDF  for  the  produet 


Mpller  state  on  the  \  ]  >  surfaee  remains  lower  for  a  wider  range  of  7  values.  On  the  other 


hand,  if  the  reaetant  Mpller  state  starts  on  the  |  ]  >  surfaee,  probability  ean  be  transferred  to 


the  1 1  >  via  radial  eoupling  at  low  values  of  7.  Although  the  effeetive  potential  eventually 


eauses  the  probability  transfer  to  the  |  J  >  to  turn  off,  the  amount  of  probability  going  to  the 


1 1  >  surfaee  is  signifieant  eompared  to  the  total  probability.  For  this  reason,  approximately 

2 

half  the  states  pumped  into  the  do  not  eontribute  to  the  relaxation  of  alkali  metals  to 
the  ^Pi_  where  lasing  oeeurs. 
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Figure  6.72:  The  maximum  value  of  the  seattering  Wigner  Distribution  Funetion  for  K+Ne 

1  3 

system  with  reaetant  M0ller  state  in  li  >,  produet  Mpller  state  in  li  >  (Blue),  and  produet 

2  2 
3 

Mpller  state  I3  >  (Red).  The  maximum  value  was  sampled  for  every  AJ  =  5  ranging  from 

2 

7  =  5.5  to  7  =  250.5  (ineluding  7  =  1.5).  Note  that  this  is  the  same  plot  shown  in  Fig.  6.49. 
This  plot  is  shown  again  as  a  matter  of  eonvenienee  for  the  reader. 


Figure  6.73:  The  maximum  value  of  the  seattering  Wigner  Distribution  Funetion  for  K+Ar 

1  3 

system  with  reaetant  Mpller  state  in  |i  >,  produet  Mpller  state  in  li  >  (Blue),  and  produet 

2  2 
3 

Mpller  state  I3  >  (Red).  The  maximum  value  was  sampled  for  every  A7  =  5  ranging  from 
7  =  5.5  to  7  =  250.5  (ineluding  7  =  1.5). 


153 


X  10 


o  o 

O  O  O  <M 

CD  <M  O  I  I 


Figure  6.74:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  1 

M0ller  state  in  |  i  >’  product  Mpller  state  in  I,  >,  and  7  =  1.5 
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Figure  6.75:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  1 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.76:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  1 

M0ller  state  in  |  i  >’  product  Mpller  state  in  |  i  >’  and  J  =  200.5.  The  scale  for  this 

2  2 

figure  is  different  than  the  scale  used  in  Fig.  6.74.  The  scale  is  based  on  the  maximum  and 
minimum  value  of  the  scattering  WDF  to  show  that  the  amplitude  exists. 
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Figure  6.77:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.78:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  | ,  >,  and  7  =  100.5 
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Figure  6.79:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  200.5 
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Figure  6.80:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  300.5 
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Figure  6.81:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  | ,  >,  and  J  =  400.5 
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Figure  6.82:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.83:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.84:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  I3  >,  and  J  =  200.5 
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Figure  6.85:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  300.5 
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Figure  6.86:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  400.5 
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Figure  6.87:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  1 

M0ller  state  in  I3  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.88:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  1 

M0ller  state  in  I3  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.89:  Scattering  Wigner  Distribution  Function  for  K+Ar  system  with  reactant 

3  1 

M0ller  state  in  |  3  >’  product  Mpller  state  in  |  i  >’  and  J  =  200.5.  The  scale  for  this 

2  2 

figure  is  different  than  the  scale  used  in  Fig.  6.87.  The  scale  is  based  on  the  maximum  and 
minimum  value  of  the  scattering  WDF  to  show  that  the  amplitude  exists. 
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6.2.5  Alkali  Metals. 


If  the  alkali  atom  is  changed  while  the  noble  gas  atom  remains  the  same,  the  angular 
momentum  matrix  in  Eq.  6.2  will  be  affected.  The  difference  in  the  various  alkali  metals 
is  the  spin-orbit  split.  As  the  alkali  metal  increases  in  mass,  the  spin-orbit  split  grows 
larger.  A  larger  spin-orbit  spliting  means  that  higher  energies  are  needed  to  achieve  the 

1  3 

transition  from  the  \\  >  to  the  |i  >  surfaces.  In  effect,  the  probability  of  making  the 

2  2 

1  3 

transition  decreases.  This  implies  that  the  transfer  of  probability  between  the  |  \  >  and  I?  > 

2  2 

surfaces  decreases  since  the  strength  of  the  radial  coupling  decreases.  Figs.  6.31-6.48  and 
Figs.  6.90-6.120  show  the  scattering  WDF  for  the  K  Ne,  Rb  -l-  Ne,  and  Cs  +  Ne  systems 
over  a  range  of  AJ  =  50  from  J  =  50.5  to  7  =  250.5  (7  =  1.5  is  included)  with  the  reactant 

3 

Mpller  state  starting  on  the  \]  >  surface.  For  the  K  +  Ne  system,  the  spin-orbit  split  is 

2 

relatively  small,  so  the  radial  coupling  is  relatively  strong.  As  such  we  see  that  there  is 

significant  probability  on  all  three  surfaces  for  the  certain  values  of  7.  The  Rb  +  Ne  system 

1 

has  a  larger  spin-orbit  split,  so  the  probability  of  making  a  transition  to  the  |  ]  >  exists  but 

2 

1 

smaller  than  the  K  +  Ne  case.  The  amplitude  of  the  scattering  WDF  for  K  +  Ne  for  the  |  ?  > 

2 

,  3 

surface  at  7  =  1.5  is  approximately  of  the  amplitude  of  the  scattering  WDF  for  the  I?  > 

2 

surface  at  7  =  1.5.  In  comparison,  the  change  in  the  order  of  magnitude  between  these  two 
surfaces  for  Rb  +  Ne  at  J  =  1.5  is  joo^ooo-  spin-orbit  split  for  Cs  is  much  larger,  so  the 
radial  coupling  term  is  very  weak.  For  this  reason,  the  change  in  the  order  of  magnitude 
in  the  scattering  WDF  between  these  surfaces  for  Cs  -l-  Ne  is  10“  10.  With  such  a  drastic 
change  in  the  amplitude  of  the  scattering  WDF  between  the  two  surfaces  for  C5  -l-  Ne,  we 
can  assume  that  there  is  almost  no  radial  Coupling  for  Cs.  For  the  purpose  of  the  DPAF, 
this  is  detrimental  to  optimizing  the  output  power  of  the  laser;  however,  for  purposes  of 
studying  a  two  level  system  that  only  explores  Coriolis  coupling,  this  can  be  beneficial. 
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Figure  6.90:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  1 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.91:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  1 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  50.5 
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Figure  6.92:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  1 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.93:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  1 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  150.5 
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Figure  6.94:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.95:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  50.5 
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Figure  6.96:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.97:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 
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M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  150.5 
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Figure  6.98:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 
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M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  200.5 
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Figure  6.99:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  250.5 
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Figure  6.100:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.101:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 

3  3 

M0ller  state  in  |  j  >,  product  Mpller  state  in  ||  >,  and  J  =  50.5 
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Figure  6.102:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 
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M0ller  state  in  |  j  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.103:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 
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M0ller  state  in  |  j  >,  product  Mpller  state  in  ||  >,  and  J  =  150.5 
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Figure  6.104:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 
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M0ller  state  in  |  j  >,  product  Mpller  state  in  ||  >,  and  J  =  200.5 
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Figure  6.105:  Scattering  Wigner  Distribution  Function  for  Rb+Ne  system  with  reactant 
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M0ller  state  in  |  j  >,  product  Mpller  state  in  ||  >,  and  J  =  250.5 
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Figure  6.106:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  1 

M0ller  state  in  |  i  >’  product  Mpller  state  in  I,  >,  and  J  =  1.5.  Note  that  most  of  the  detail 

2  2 

in  this  figure  is  numerical  error  associated  with  the  Fourier  transform.  The  error  is  visible 
because  the  scale  is  on  the  order  of  10“^. 
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Figure  6.107:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  I 

M0ller  state  in  |  i  >’  product  Mpller  state  in  | ,  >,  and  J  =  50.5.  Note  that  most  of  the  detail 

2  2 

in  this  figure  is  numerical  error  associated  with  the  Fourier  transform.  The  error  is  visible 
because  the  scale  is  on  the  order  of  10“^. 
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Figure  6.108:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  1 

M0ller  state  in  |i  >,  product  Mpller  state  in  I,  >,  and  J  =  100.5.  Note  that  most  of  the 

2  2 

detail  in  this  figure  is  numerical  error  associated  with  the  Fourier  transform.  The  error  is 
visible  because  the  scale  is  on  the  order  of  10“^. 
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Figure  6.109:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 
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M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.110:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 
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M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  50.5 
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Figure  6.111:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 
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M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.112:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  150.5 
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Figure  6.113:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  200.5 
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Figure  6.114:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  250.5 
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Figure  6.115:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  1.5 
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Figure  6.116:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  50.5 


2 


2 


197 


Time  (A.U) 


X  10 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

C\J 

CO 

00 

■!— 

00 

CO 

C\J 

o 

1 

1 

1 

1 

CO 

I 


un 

O) 


00 


CD 


in 


CO 


cj)oor^cDLn'^coc\j-i-o 

XSiaug 


T  I  I  r 


J _ I _ I _ L 


Figure  6.117:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  7  =  100.5 
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Figure  6.118:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  150.5 
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Figure  6.119:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 

3  3 

M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  200.5 
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Figure  6.120:  Scattering  Wigner  Distribution  Function  for  Cs+Ne  system  with  reactant 
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M0ller  state  in  |i  >,  product  Mpller  state  in  ||  >,  and  J  =  250.5 
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6.2.6  Negative  Chirp. 

Like  the  square  well  problem,  we  are  interested  in  the  shape  of  the  negative  Chirp 
that  the  seattering  WDF  makes  in  the  time-energy  axis.  We  to  get  a  rough  idea  of  the 
elassieal  trajeetory  that  a  partiele  with  a  eertain  energy  makes  while  under  the  influenee 
of  the  Hamiltonian  defined  in  Eq.  6.2.  From  the  adiabatie  potential  energy  surfaees,  we 
assume  that  we  have  a  partiele  in  a  zero-potential  field.  This  partiele  will  eollide  elastieally 
with  an  infinitely  high  wall.  The  loeation  of  the  wall  is  determined  by  the  interseetion  of 
the  peak  momentum  of  the  reaetant  Mpller  state  with  the  adiabatie  potential  energy  surfaee. 

The  potential  energy  surfaee  used  is  the  surfaee  at  whieh  the  reaetant  Mpller  state  starts  on. 

1 

For  this  example,  we  will  use  the  |  ]  >  adiabatie  surfaee.  Fet  d  be  the  distanee  between  the 

2 

loeation  of  the  wall  and  the  loeation  of  reaetant  Mpller  state.  Beeause  the  elassieal  partiele 
experienees  no  potential,  the  total  energy  is  the  kinetie  energy  whieh  is  F  =  For  a 

given  veloeity  v,  the  time  the  partiele  takes  to  travel  to  the  wall  and  baek  to  the  starting 
loeation  (whieh  is  where  the  produet  Mpller  state  is  loeated)  is  given  by  t  In  this 

model,  the  energy  time  relation  is  given  by 


Although  this  model  is  simple,  the  square  well  and  barrier  like  features  on  the  adiabatie 
PES  will  have  a  perturbing  effeet  on  Eq.  6.4.  Using  the  same  Matlab  imaging  statistieal 
analysis  teehniques  that  were  used  for  the  square  well  ehirp,  we  ean  arrive  at  an  energy¬ 
time  relation.  A  power  eurve  of  the  form  E  =  at^,  where  a,  b,  and  c  are  fit  parameters, 
provided  the  best  fit  with  b  =  -2.119  +  0.001  on  a  99  %  eonfidenee  interval.  The  same 

13  3  3 

eurve  was  also  plaeed  on  top  of  the  seattering  WDF  assoeiated  with  cX\(t)  and  cXlif) 

2’2  2’2 

eorrelation  funetions.  The  values  of  the  seattering  WDF  for  eaeh  eorrelation  funetion  were 
plotted  on  the  same  seale  in  order  to  see  the  overlap  of  peaks  within  eaeh  band  as  shown 
in  Fig.  6.121.  This  ealeulation  was  repeated  for  multiple  bands.  At  high  energies  and  short 
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times,  the  peaks  are  relatively  out  of  phase  at  low  energies  and  long  times,  the  peaks  are 
relatively  in  phase. 


Figure  6.121:  Amplitude  of  the  Seattering  WDF  for  the  produet  Mpller  states  on  the  I?  > 

2 

3 

(blue  line)  and  the  I3  >  (red  line)  surfaees  along  the  ehirp  at  J  =  100.5  for  the  K  +  Ne 

2 

]_ 

system.  The  shape  of  the  ehirp  was  determined  from  the  seattering  WDF  for  the  |  ?  > 

2 

surface.  The  x-axis  corresponds  to  a  point  along  the  chirp  while  the  y-axis  is  the  amplitude 
of  the  scattering  WDF. 


203 


VII.  Conclusion 


The  scattering  Wigner  Distribution  Function  was  able  to  give  us  some  insight  into  the 
alkali  metal  plus  noble  gas  collision.  The  square  well  scattering  WDF  provided  insight  into 
interpreting  the  meaning  behind  the  positive  and  negative  values  of  the  scattering  WDF. 
Due  to  the  property  of  the  scattering  WDF  that  the  integrated  projection  must  be  the  related 
to  the  s-matrix  element  up  to  a  normalization  factor,  the  positive  values  corresponds  to  the 
presence  of  probability  at  a  particular  energy  in  a  moment  in  time.  The  negative  values 
correspond  to  the  loss  of  probability  at  that  surface  at  a  particular  energy  in  a  moment 
of  time.  The  location  of  this  ’’lost”  probability  most  likely  ended  up  in  the  reflection 
coefficient,  but  further  analysis  must  be  done  before  this  can  be  confirmed.  The  zero  value 
of  the  scattering  WDF  means  that  the  probability  was  elsewhere  in  the  system. 

Using  this  interpretation  for  the  meaning  behind  the  values  of  the  scattering  WDF,  we 
were  able  to  comprehend  the  meaning  behind  the  scattering  WDF  for  the  M  +  Ng  collision. 
As  the  total  angular  momentum  J  increased,  the  Coriolis  coupling  and  effective  potential 
grew  stronger  while  the  strength  of  the  radial  coupling  remained  constant.  If  the  reactant 

3 

Mpller  state  started  on  the  \]  >  surface,  there  was  a  value  of  J  such  that  the  transfer  of 

2 

3 

probability  to  the  li  >  surface  was  cut  drastically  reduced  compared  to  when  J  =  1.5. 

2 

As  the  mass  of  the  noble  gas  partner  increased,  the  value  of  J  at  which  the  probability  on 
1 

the  Ij  >  surface  began  to  diminish  increased.  Although  the  peak  amount  of  probability 

2 

1 

transferred  to  the  |  j  >  was  diminished  for  heavier  noble  gas  atoms  than  lighter  atoms,  the 

2 

range  of  total  angular  momentum  for  which  there  were  significant  probability  transfer  was 
much  greater  for  the  heavier  noble  gas  atoms. 

Future  work  can  focus  on  further  investigating  the  correlation  between  the  number  of 
bands  with  oscillating  peaks  and  valleys  and  the  amplitude  of  the  scattering  WDF.  From  the 
square  well,  we  see  that  the  number  of  bands  with  oscillating  peaks  and  valleys  increased 
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as  the  well  depth  increased  which  corresponded  to  an  increase  in  probability  of  being 
transmitted  through  the  square  well  on  the  secondary  transmissions.  The  appearance  of 
these  bands  of  oscillating  peaks  and  valleys  in  the  M  +  Ng  system  correlated  with  the 
amplitude  of  the  scattering  WDF.  The  dip  like  structure  that  was  found  in  the  amplitude  of 

3 

the  scattering  WDF  for  the  product  Mpller  state  on  the  1 1  >  surface  corresponded  to  the 

2 

1 

absence  of  these  bands.  The  amplitude  of  the  scattering  WDF  for  the  \  \  >  surface  should 

2 

be  plotted  in  order  to  see  if  there  is  a  corresponding  peak  at  this  surface.  Work  was  focused 
on  the  frequency  at  which  the  bands  present  in  the  scattering  WDF  for  the  M  +  Ng  arrived 
at  the  product  Mpller  state;  however,  perhaps  the  frequency  of  bands  with  oscillating  peaks 
and  valleys  or  the  number  of  oscillating  peaks  and  valleys  should  be  focused  on  instead. 
Future  work  could  also  focus  on  different  phase  space  representations  of  the  scattering 
event.  The  Wigner  Distribution  Function  belongs  to  the  Cohen  class  of  distributions [2,  15]. 
One  such  function  that  can  be  studied  is  the  ambiguity  function.  The  ambiguity  function 
is  the  double  Fourier  transform  of  the  Wigner  Distribution  Function  with  respect  to  delay 
time  and  energy  values,  r  and  Future  work  can  involve  applying  the  scattering  WDF  to 
various  other  collisional  systems  of  interest  or  could  involve  more  detailed  studies  of  the 
current  systems  of  interest.  For  example,  the  scattering  WDF  could  be  computed  for  the 
coupling  between  higher  energy  states  in  the  H  +  H2,  OH  +  CO,  and  M  +  Ng  systems. 
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